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noise  measurements,  and  provides  evidence  that  our  model  for  the  Instrument's 
response  to  gradient  fields  Is  correct./'  However,  experiments  designed  to 
Investigate  balancing  features  of  our  gradlometer  show  that  the  Instrument 
response  under  rotation  In  a  uniform  magnetic  field  Is  In  marked  disagreement 
with  the  model.  Finally  we  describe  the  two  types  of  anamolous  behavior 
observed  during  the  balancing  experiments. 


SUMMARY 


Measurements  of  fluctuating  magnetic  gradients  above  the  ocean's  surface 
provide  knowledge  of  water  motions  below  the  surface.  Using  a  superconduct¬ 
ing  magnetic  gradiometer  we  have  recently  completed  measurements  of  very  small 
fluctuating  magnetic  gradients  associated  with  internal  ocean  waves  passing  the 
Naval  Ocean  Systems  Center  (NOSC)  Oceanographic  Research  Tower  located  near 
San  Oiego,  CA.  The  extreme  sensitivity  of  the  gradiometer  and  small  ampli¬ 
tude,  low  frequency  nature  of  the  signal  of  interest  dictated  the  develop¬ 
ment  of  effective  noise  suppression  techniques.  In  this  report  we  present 
results  of  extensive  efforts  to  develop  techniques  which  would  allow  measure¬ 
ments  of  fluctuating  magnetic  gradients  of  order  1  pi  cotesla/meter  at  freq¬ 
uencies  of  10~2  to  10"3  Hz  in  the  presence  of  a  gradient  of  200  nanotesla/meter. 

First  we  present  a  chronological  review  of  experimental  work  to  date  and 
a  description  of  our  superconducting  magnetic  gradiometer.  In  Section  III 
we  describe  our  theoretical  model  for  the  gradiometer  response  in  an  arbitrary 
magnetic  field  as  the  sum  of  two  terms;  a  gradient  term  and  an  imbalance  term. 
The  imbalance  term  arises  from  small  differences  in  area  and  orientation  of 
the  gradiometer  sensing  loops  and  produces  a  magnetometer-like  response  during 
gradiometer  rotations;  the  gradient  term  gives  the  gradiometer  output  as  a 
function  of  gradiometer  orientation  and  ambient  gradient  field.  We  represent 
a  steady  gradient  field  as  arising  from  a  fixed  equivalent  dipole  (Section  IV) 
and  small  gradient  fluctuations  are  treated  as  perturbations  in  the  position 
and  orientation  of  the  equivalent  uipole.  When  combined  with  our  description 
of  the  gradiometer  response,  the  equivalent  dipole  representation  provides 
a  model  for  external  gradient  noise  produced  in  the  gradiometer  by  large  steady 
gradients  from  the  steel  structure  of  the  NOSC  tower. 


In  Section  V  we  describe  various  sources  of  gradioineter  noise  and  general 
techniques  of  noise  suppression.  We  develop  the  concept  of  a  minimum  noise 
orientation  (MNO)  in  Section  VI  in  the  context  of  the  equivalent  dipole  model 
and  our  representation  of  the  gradiometer  response  to  an  ambient  gradient  field. 
We  also  describe  an  extremely  useful  technique  for  determining  the  MNO  and  show 
results  of  measurements  in  both  a  simulated  gradient  environment  and  at  the 
NOSC  tower.  Low  frequency  gradiometer  power  spectra  demonstrate  the  effect¬ 
iveness  with  which  the  alignment  technique  suppresses  gradient  noise  in  the 
instrument. 

Section  VII  deals  with  general  gradiometer  balancing  techniques  and  des¬ 
cribes  unusual  and  unexpected  behavior  encountered  during  balancing  experiments 
with  our  present  instrument.  Large  hysteretic  effects  appear  as  the  gradiom¬ 
eter  is  tilted  more  than  30°  from  vertical.  Smaller  non-hysteretic  effects 
limit  the  accuracy  with  which  we  can  define  the  gradiometer  balance  in  the  con¬ 
text  of  our  mathematical  model  for  the  instrument's  response.  These  anamolies 
impact  significantly  on  our  ability  to  precisely  measure  an  ambient  gradient 
field  and  balance  the  gradiometer  in  the  presence  of  a  steady  gradient  (Sec¬ 
tion  VIII).  We  also  describe  various  attempts  to  obtain  precise  gradiometer 
alignment  for  noise  suppression  within  the  constraints  imposed  by  these  anamol¬ 
ies.  The  failure  of  these  techniques  eventually  led  to  development  of  our 
empirical  MNO  measurements.  In  Section  IX  we  present  results  of  a  computer 
calculation  which  accurately  describes  the  empirical  observations.  The  success 
of  the  computer  model  confirms  that  we  can  represent  small  gradient  fluctua¬ 
tions  at  the  gradiometer  as  arising  from  perturbations  in  the  position  and 
orientation  of  an  equivalent  dipole.  The  general  nature  of  the  model  suggests 
that  the  technique  may  be  applicable  to  a  variety  of  situations  involving 
magnetic  gradient  measurements  in  the  presence  of  large  ambient  fields. 
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I.  INTRODUCTION 


Seawater  moving  across  the  steady  magnetic  field  of  the  earth  generates 
electric  currents  which  give  rise  to  small  magnetic  fields  and  magnetic  grad¬ 
ients  both  below  and  above  the  ocean's  surface.  From  measurements  of  the  mag¬ 
netic  fields  or  gradients  above  the  surface  we  deduce  the  motions  of  seawater 
below  the  surface.  Commercially  available  superconducting  magnetometers  and 
gradiometers,  developed  in  recent  years  with  their  unprecedented  sensitivity 
and  low  inherent  noise,  provide  a  tool  with  which  to  make  such  measurements. 

Surface  waves,  which  produce  velocities  of  order  100  cm/sec,  generate 
magnetic  fields  and  gradients  which  are  readily  detectable  wi th  a  variety  of 
instruments.  Water  velocities  from  internal  waves  however  are  nearly  two  or¬ 
ders  of  magnitude  smaller  and  the  superconducting  sensor  provides  a  unique 
device  with  which  to  detect  magnetic  effects  from  internal  waves.  In  a  recent 
experiment  we  measured  the  fluctuating  magnetic  gradients  above  the  water 
produced  by  internal  ocean  waves  passing  the  NOSC  Oceanographic  Research  Tower. 
These  measurements,  when  combined  with  simultaneous,  direct  measurements  of 
water  displacement  and  velocity  allow  us  to  determine  an  empirical  relation¬ 
ship  between  subsurface  water  motions  and  the  associated  magnetic  effects,  a 
necessary  first  step  in  attempting  to  deduce  seawater  motions  through  measure¬ 
ments  of  magnetic  gradients  above  the  surface. 

The  primary  obstacle  to  performing  such  measurements  is  the  magnetic 
noise  introduced  into  the  gradiometer  from  effects  not  associated  with 
internal  waves.  Although  the  magnetic  gradient  spectrum  we  initially 
calculated  for  linear  internal  waves  lies  well  above  the  expected  inherent 
noise  spectrum  of  the  SQUID  gradiometer  (Podney  and  Gillespie,  1977)  to 


insure  a  high  probability  of  detecting  and  measuring  the  gradients  of  in¬ 
terest  all  other  magnetic  noise  must  be  suppressed  at  least  to  a  level 
comparable  to  the  inherent  sensor  noise.  This  problem  was  greatly  intens¬ 
ified  by  the  necessary  proximity  of  the  steel  structure  of  the  NOSC  tower, 
which  produced  large  gradients  at  the  site  of  the  gradiometer  making  the 
instrument  highly  sensitive  to  small  relative  motions  between  the  gradio¬ 
meter  and  tower.  In  addition  the  gradients  generated  by  the  tower  fluctuate 
in  response  to  changes  in  the  earth's  field  (as  a  result  of  ionospheric 
activity)  producing  magnetic  noise  at  internal  wave  frequencies  even  in 
an  absolutely  rigid  tower-gradiometer  system.  Other  less  serious  con¬ 
tributions  to  the  magnetic  environment  of  the  gradiometer 
further  degrade  the  effective  signal-to-noise  ratio.  Hence,  prior  to 
attempting  to  measure  the  magnetic  effects  from  internal  waves  it  was 
necessary  to  investigate  the  capabilities  and  limitations  of  our  super¬ 
conducting  gradiometer  and,  through  field  testing,  develop  techniques  to 
provide  sufficient  suppression  of  the  noise  from  nearby  magnetic  objects 
to  allow  the  measurements  to  be  performed  from  the  NOSC  tower. 

This  document  reports  the  results  of  experiments  at  the  NOSC  La  Posta 

*  ** 
Astrogeophysical  Observatory  and  at  the  NOSC  Oceanographic  Research  Tower 

*The  La  Posta  Astrogeophysical  Observatory  is  located  in  the  mountains  approx¬ 
imately  70  miles  east  of  San  Diego,  California  at  116°  25*  6"  west  longitude 
and  32°  40'  39"  north  geodetic  latitude  (Bleiweiss  and  Hefer,  1975). 

The  tower  is  located  approximately  one  mile  off  Mission  Beach  near  San  Diego, 

California  in  about  18  meters  of  water. 
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designed  to  develop  the  requisite  noise  suppression  techniques  and  investigate 
the  inherent  properties  of  our  instrument.  The  mathematical  basis  for  most 
of  the  noise  suppression  experiments  described  herein  was  presented  in  a  prev¬ 
ious  report  as  were  initial  measurements  of  instrument  noise  and  instrument 
balance  (Podney  and  Gillespie,  1977).  For  most  of  this  information  we  refer 


to  that  earlier  report  and  include  here  only  the  brief  descriptions  required 
for  a  clear  understanding  of  this  presentation. 


II.  EXPERIMENTAL  BACKGROUND 


To  place  the  work  described  In  this  report  into  its  proper  perspec¬ 
tive  from  the  viewpoint  of  the  overall  experimental  objectives  we  present 
in  this  section  a  brief  outline  of  the  experimental  development  to  date 
including  a  description  of  the  Instrument  and  modifications  made  in  March, 
1978. 

A.  Chronological  Review 

Field  tests  with  the  instrument  began  in  January  of  1976  shortly 
after  wo  acquired  the  instrument,  and  preliminary  measurements  of  instru¬ 
ment  noise  spectra,  balancing  capabilities,  and  gradiometer  alignment  in 
a  steady  gradient  field  continued  throughout  19/6.  Those  results  were 
reported  in  Podney  and  Gillespie,  1977.  Early  balancing  and  alignment 
experiments  which  used  a  wooden  gimbal  to  rotate  and  orient  the  gradiom¬ 
eter,  were  not  conclusive  as  the  gimbal  was  limited  to  only  two  indepen¬ 
dent  axes  of  rotation.  As  discussed  in  detail  In  Podney  and  Gillespie, 

1977,  our  proposed  alignment  and  noise  suppression  techniques  required 
the  capability  to  rotate  the  Instrument  about  three  Independent  axes  of 
rotation.  Consequently,  further  tests  were  delayed  until  the  acquisition 
of  our  present,  fiberglass,  3-axis,  gimbal  on  July  ?S,  1977.  Figure  la  is 
a  sketch  of  the  gimbal.  From  this  date  until  about  August  16,  1977  further 
field  tests  were  conducted  at  the  la  Posta  site  to  complete  the  development 
of  the  techniques  to  balance  and  align  the  gradiometer  in  a  large  gradient 
field  and  suppress  the  noise  from  nearby  magnetic  objects.  At  the  completion 
of  these  tests  the  gradiometer  was  transported  to  the  NOSC  oceanographic  tower 


ROTATION  ABOUT 
CRYOSTAT  AXIS 


ROTATION  ABOUT 
VERTICAL  AXIS 


GIBMAl 


TILING  ABOUT 
HORIZONTAL  AXIS 


Figure  la.  Sketch  of  gimbal  mount  holding  cryostat  at  an  angle  of  30°  from 
vertical.  Tilting  about  a  horizontal  axis  (limited  to  +45°)  and 
two  axes  of  rotation  (one  vertical  and  the  other  directed  along 
the  axis  of  the  cryostat)  allow  every  orientation  of  pickup  loops 
forming  the  gradiometer  within  a  90°  solid  angle  about  vertical. 


where  internal  wave  measurements  were  performed  from  mid-August 
through  the  first  week  of  October.  The  results  of  these  measurements 
(which  will  be  discussed  in  a  separate  report)  were  inconclusive  in 
themselves  but  suggested  that  further  measurements  with  only  slightly 
improved  noise  suppression  could  produce  excellent  quality  data  with 
which  mathematical  predictions  could  be  compared. 

In  November,  1977  we  resumed  field  tests  at  the  La  Posta  site 
which  continued  through  the  end  of  April  1978  (with  the  exception  of 
March  when  the  instrument  was  warmed  up  and  several  design  modifica¬ 
tions  were  introduced).  During  this  period  we  developed  and  tested 
an  empirical  technique  for  aligning  the  gradiometer  in  the  steady 
gradient  field  of  an  iron  sphere  which  was  used  to  simulate  the  gradi¬ 
ent  field  from  the  steel  tower  and  demonstrated  that  with  proper  align¬ 
ment,  the  magnetic  noise  from  the  sphere  could  be  suppressed  to  below 
the  level  of  inherent  instrument  noise  over  the  frequency  band  of  interest. 
Although  less  successful  in  our  attempts  to  balance  the  gradiometer 
in  the  presence  of  a  large  gradient,  we  did  develop  a  mechanical  means 
to  insure  a  level  of  balance  adequate  to  suppress  the  noise  introduced 
directly  into  the  gradiometer  by  the  uniform  earth's  field. 

In  May,  1978  we  again  transported  the  gradiometer  to  the  NOSC 
tower  where  we  demonstrated  the  new  noise  suppression  techniques  in 
an  oceanographic  environment  but  in  the  relative  absence  of  internal 
wave  activity.  The  success  of  these  trials,  which  continued  through 
the  first  week  of  June,  provided  the  motivation  for  a  second  attempt 
at  the  internal  wave  measurements  which  subsequently  were  conducted 


from  17  July  through  24  August  1978.  The  results  of  the  1978  internal  wave 
measurements  will  also  be  presented  in  the  separate  report  referred  to  above. 

The  present  report  is  concerned  primarily  with  the  work  performed  between 
November  1,  1977  and  June  30,  1978,  summarized  by  the  following  tasks: 

a.  Made  several  design  modifications  to  the  gradiometer  probe  assembly 
during  March  1978,  and  readjusted  the  permanently  affixed  balance  disks 
following  the  changes. 

b.  Remeasured  instrument  noise  spectra  after  making  instrument  modif¬ 
ications. 

c.  Developed  and  tested  procedures  for  aligning  the  gradiometer  in  a 
steady  gradient  field  to  maximize  noise  suppression  from  nearby  mag¬ 
netic  objects  (using  an  iron  sphere  at  La  Posta). 

d.  Attempted  to  develop  a  technique  to  balance  the  instrument  in  the 
presence  of  a  large  gradient  field.  Although  not  successful  in  this 
endeavor,  we  did  develop  means  to  insure  an  adequate  level  of  bal¬ 
ance  on  the  tower. 

e.  Demonstrated  the  noise  suppression  alignment  technique  at  the  NOSC 
tower. 

In  the  following  sections  we  will  describe  in  detail  the  results  of  this  work 
but  first  we  present  a  brief  description  of  the  instrument  and  the  modifica¬ 
tions  we  introduced  to  help  improve  certain  aspects  of  its  performance. 

B.  Instrument  Description 

*  .  . 

Our  gradiometer,  which  was  purchased  commercially,  is  a  single  axis 
device  containing  two  coplanar  pick-up  loops  with  their  centers  ap¬ 
proximately  25  cm  apart.  The  counterwound  loops,  each  4.45  cm  in 

*The  instrument  was  built  by  Superconducting  Technology,  Inc.  (SCT). 


diameter,  are  constructed  of  superconducti ng  wire  and  are  connected 
to  a  Superconducting  Quantum  Interference  Device  (SQUID)  via  a  Radio 
Frequency  Interference  (RFI)  isolation  transformer.  The  loops  are 
glued  into  circular  grooves  machined  in  a  high  purity  silicon  substrate 
which  provides  rigidity,  a  mechanical  base  for  mounting,  and  helps 
reduce  thermal  gradients  in  the  helium  bath.  The  entire  assembly  of 
SQUID,  substrate,  and  sensing  loops  is  then  supported  by  a  probe  assembly 
which  is,  in  turn,  inserted  into  a  socket  in  the  bottom  of  the  cryogenic 
dewar.  The  socket  serves  to  position  and  hold  the  probe  assembly  in 
place.  The  probe  assembly  accomodates  three  pushrods  with  which  to  adjust 
the  small  niobium  disks  used  to  balance  the  instrument,  and  also  has 
provision  for  the  radio  frequency  (RF)  transmission  line  via  which  the 
RF  and  audio  bias  signals  are  applied  to  the  SQUID  and  the  SQUID  output 
signal  is  received.  Figure  lb  shows  the  overall  assembly  including  the  cry¬ 
ogenic  dewar.  The  top  of  the  probe  assembly  provides  for  mounting  the 
RF  electronics  used  to  bias  the  SQUID  sensor  and  a  three  axis 
fluxgate  magnetometer.  A  more  detailed  description  of  the  instrument 
and  its  design  features  can  be  found  in  Podney  and  Gillespie,  1977. 

C.  Instrument  Modifications 

At  least  part  of  the  success  of  the  May  and  July,  1978  experiments 
at  the  NOSC  tower  must  be  attributed  to  the  improvement  in  the  inherent 
noise  levels  of  the  gradiometer  system  following  the  instrument  modifi¬ 
cations  introduced  during  March  1978.  During  balancing  trials  in  Jan¬ 
uary  1978  it  was  extremely  difficult  to  adjust  the  balance  disks,  partly 
as  a  result  of  solidified  nitrogen,  oxygen,  and  carbon  dioxide  building 
up  on  the  probe  assembly  inside  the  dewar.  Furthermore,  we  observed 
sudden  large  dc  shifts  in  the  gradiometer  output  as  the  instrument  was 
tilted  to  angles  greater  than  about  30°,  suggesting  some  type  of  mechan- 
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Figure  lb.  Illustration  of  the  gradiometer  and  cryogenic  container  showing 

the  principle  conponents  of  the  gradiometer  probe.  A  wrapping  of 
alternating  layers  of  fiberglass  cloth  and  aluminized  mylar  (not 
shown)  insulates  the  interior  vessel, 
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leal  motion  occuring  inside  the  dewar.  With  this  motivation  we  trans¬ 
ported  the  gradiometer  and  dewar  to  temporary  laboratory  facilities  at 
NOSC  and  removed  the  gradiometer  from  the  dewar  after  warming  the  entire 
system  up  to  room  temperature.  Our  initial  examination  of  the  instrument 
revealed  the  following: 

a.  The  difficult  operation  of  the  balance  slides  apparently  resulted 
from  two  factors.  The  0-ring  seal  at  the  top  plate  of  the  probe 
had  become  dry  and  could  not  be  lubricated  without  removing  the 
probe  from  the  dewar  and  the  feed-through  design  at  the  top  plate 
allowed  dirt  to  accumulate  above  the  0-ring  seal  and  cause  binding 
of  the  balance  rods.  A  build  up  of  solid  material  on  the  portion 
of  the  probe  held  at  cryogenic  temperatures  undoubtedly  contrib¬ 
uted  further  to  the  problem. 

b.  The  instrument  was  held  in  place  in  the  Kevlar  support  assembly 
at  the  bottom  of  the  probe  with  four  stainless  steel  screws  all 

of  which  were  within  two  inches  of  the  gradiometer  pick  up  loops. 

As  verification  of  the  magnetic  nature  of  the  screws,  when  suspended 
with  a  thread  about  20  inches  long  the  screws  could  be  deflected 
about  3  inches  from  vertical  with  a  small  magnet. 

c.  A  small  rod  of  G-10  fiberglass  .8  mm  in  diameter  which  provided 
tension  on  the  balance  slides  was  cracked  partly  releasing  the 
tension  on  one  of  the  balance  slides. 

d.  The  battery  operated  helium  level  indicator  was  not  functioning. 

The  following  modifications  were  made  to  the  gradiometer  and  probe  assembly 
to  correct  the  problems  listed  above. 

a.  The  top  plate  feed  throughs  were  redesigned  to  allow  proper  lub¬ 
rication  of  the  0-ring  seal  without  opening  the  dewar  and  to 
eliminate  the  area  in  which  dirt  could  collect.  Also  protective 
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r. 


covers  for  the  tops  of  the  balance  rods  were  constructed  to  help 
k.  maintain  cleanliness  around  the  balance  rod  seals  during  field 

operations.  Finally,  the  lower  end  of  the  balance  rods  were 
redesigned  to  provide  a  bayonet-type  coupling  mechanism  with  which 
the  balance  rods  could  be  completely  decoupled  from  the  balance 
slides.  A  vertical  travel  limiter  and  guide  for  the  bayonet  coupling 
was  incorporated  into  the  new  feed-through  design  at  the  top  plate. 

b.  The  stainless  steel  screws  were  replaced  with  screws  made  of 
G-10  fiberglass. 

c.  The  .8  mm  G-10  fiberglass  rod  was  replaced  to  provide  positive 
tension  on  all  of  the  balance  slides. 

d.  The  battery  operated  helium  level  indicator  was  removed  and  we 
subsequently  relied  on  a  manual  dipstick  to  measure  the  helium 
level  in  the  dewar. 

The  design  changes  made  in  the  balance  rods  worked  well  and  no  further 
binding  of  the  balance  rods  was  encountered.  The  bayonet  coupling  scheme 
was  of  only  marginal  use  due  to  slight  misalignment  in  the  original  probe 
assembly  making  it  difficult  to  reengage  one  of  the  balance  slides  once 
the  rod  had  been  disengaged.  This  feature  was  not  used  beyond  some  limited 
testing  following  the  initial  cool  down  of  the  gradiometer.  The  replace¬ 
ment  of  the  stainless  steel  screws  with  nonmagnetic  screws  of  G-10  fiber¬ 
glass  presented  a  more  serious  problem.  After  returning  to  La  Posta  to 
resume  field  tests  we  immediately  found  that  we  could  no  longer  balance 
the  gradiometer.  More  specifically,  as  a  result  of  removing  the  stain¬ 
less  steel  screws,  the  gradiometer  imbalance  was  so  large  that  the  small 
adjustable  niobium  balance  disks  could  not  compensate  for  it.  Consequently 
it  was  necessary  to  modify  the  permanently  affixed  balance  disks  on  the 
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gradiometer  which  required  that  the  instrument  be  removed  from  the  dewar 
and  allowed  to  warm  up  to  room  temperature.  This  procedure  was  performed 
twice  before  we  were  able  to  restore  the  permanent  balance  of  the  instru¬ 
ment  to  a  level  which  would  allow  compensation  using  the  adjustable  bal¬ 
ance  disks. 

Before  making  the  modifications  described  above  we  performed  a  measure¬ 
ment  of  the  inherent  instrument  noise  spectrum  at  the  La  Posta  site. 

After  the  modifications  had  been  completed  we  returned  the  instrument  to 
La  Posta  and  repeated  the  measurement.  The  power  spectra  derived  from 
these  two  measurements  are  shown  in  Figure  2.  The  heavy  solid  line  shows 
the  SQUID  sensor  noise  measured  previously  (Podney  and  Gillespie,  1977), 
the  broken  line  shows  the  instrument  noise  measured  before  the  modifications 
and  the  dashed  line  shows  the  instrument  noise  after  the  modifications  were 
made.  In  the  frequency  range  of  1  to  5  milliHertz  (mHz)  the  instrument 
noise  power  was  reduced  by  about  a  factor  of  3,  which  can  probably  be 
accounted  for  by  removal  of  the  stainless  steel  screws  and  any  solidified 
air,  particularly  oxygen  which  displays  temperature  dependent  paramagnetic 
behavior  at  liquid  helium  temperatures.  After  achieving  the  instrument 
noise  level  shown  in  Figure  2  which  was  comparable  to  that  originally  ex¬ 
pected  from  the  gradiometer,  we  proceeded  to  develop  techniques  that  would 
allow  us  to  obtain  similar  noise  levels  even  when  operating  the  gradiometer 
close  to  large  magnetic  objects.  This  work  is  discussed  in  the  following 


III.  THEORETICAL  MODEL  FOR  INSTRUMENT  RESPONSE 


Before  describing  the  experimental  work  concerned  with  balancing  the 
instrument  and  the  field  tests  in  which  we  developed  alignment  and  noise 
suppression  techniques  we  first  need  to  establish  the  mathematical  model 
upon  which  this  work  was  based.  Detailed  discussions  of  virtually  every 
aspect  of  this  model  were  presented  in  Podney  and  Gillespie,  1977  and  the 
detail  will  not  be  reproduced  nere.  However,  a  firm  grasp  of  the  physical  prin¬ 
ciples  behind  that  work  is  essential  to  understanding  later  sections  of  this 
report  and  we  will  refer  extensively  to  the  results  presented  there.  To 
avoid  confusion  we  will  comply  with  their  notation  as  nearly  as  possible  with 
the  following  exception. 

Throughout  the  rest  of  this  report  we  will  generally  wish  to  specify 
the  elements  of  the  gradient  matrix,  G,  with  respect  to  a  coordinate  system 
fixed  with  respect  to  the  earth.  Consequently  we  will  adopt  the  following 
notational  convention.  Unless  explicitly  stated  otherwise,  the  gradient 
matrix  G  (and  its  elements,  g.  .)  when  used  without  prime  or  superscript  will 

•  J 

represent  gradients  specified  with  respect  to  the  earth's  basis,  and  when 
primed  (G1  and  g!  )  they  will  represent  gradients  specified  with  respect 

*  J 

to  the  gradiometer  basis.  Notation  using  other  subscripts  or  superscripts 
will  be  defined  as  it  is  introduced. 

A.  Gradiometer  Response  in  a  Magnetic  Field 

To  describe  the  gradiometer  response  to  some  arbitrary  magnetic 
field  b(r,t)  which  is  a  function  of  position,  r,  and  time,  t,  we  must 
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calculate  the  differential  flux  threading  the  coplanar  counterwound  pick¬ 
up  loops.  In  some  arbitrary  coordinate  system  define  the  midpoint  of  the 

loops  to  be  at  the  position  r,  the  centers  of  the  two  loops  at  r  +  s  and 

->  ->  A  A 

r  -  s,  and  the  normals  to  the  loops  as  n^  and  n^.  The  flux  through  each 
loop  is  then  given  by: 


*1  =  /  n^b  (r+s+pj  ,t) 

A1 


0) 


$2  =  J  dp2  n^-t)  (r-s+p2,t) 
A2 


(2) 


where  p-j  and  p2  locate  points  within  loop  contours  and  integrals  extend  over 
the  area  of  each  loop.  As  shown  in  Podney  and  Gillespie,  1977  if  the  expres¬ 
sions  are  expanded  in  a  Taylor  series  about  the  midpoint  of  the  loops,  r,  to 
first  order  the  differential  flux,  $  =  <?-j  -  ^  is  given  by 


=  2sA2  [(s-V)  n2-b  ( r , t ) ]  +  A2  [^-  n-, ~n2]  (r,t) 


(3) 


For  an  imperfectly  balanced  gradiometer  in  which  the  areas  and  orienta¬ 
tions  of  the  loops  are  slightly  different,  we  use  loop  2  as  the  reference 

A  A  A  A 

and  take  n2  =  n,  A2  =  A  to  get  6  = (l/2s  )[A^/A2  n^-n].  The  gradiometer 

A  A 

response  r(s,n,^-,t)  can  then  be  written 


r(s,n,  r,t) 


jgr*  (it)n-S(r.t)  * 


(4) 
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where  the  first  term  gives  the  response  of  the  gradiometer  to  the  grad¬ 
ient  of  the  field  at  the  midpoint  of  the  loops  while  the  second  term, 
is  proportonal  to  the  magnetic  field  at  the  midpoint.  Hence,  in 
this  picture,  the  imbalance  term  has  the  same  form  as  the  response  of 
a  vector  magnetometer  located  at  the  midpoint  of  the  loops. 

Inherent  in  the  above  discussion  is  the  assumption  that  the  grad¬ 
ients  are  essentially  constant  over  the  distance  separating  the  pick-up 
loops  so  that  the  instrument  behaves  as  a  point  gradiometer.  On  the  basis 
of  equation  (4)  the  gradiometer  output  is  treated  as  the  sum  of  two  com¬ 
ponents,  one  of  which  arises  from  the  gradient  field  as  detected  by  the  grad¬ 
iometer  and  the  other  from  the  uniform  component  of  the  field  as  a  result 
of  the  gradiometer  imbalance.  In  this  model,  once  the  gradiometer  is 
properly  balanced,  say  to  1  part  in  106,  an^  rotation  of  the  device  in 
the  earth's  field  (in  a  magnetically  clean  environment)  should  produce  a 
maximum  signal  of  order  200  pT/m  which  arises  from  the  imbalance.  (The 
gradients  in  the  earth's  field  due  to  its  dipole  nature  are  expected  to 
be  of  order  10  pT/m.)  This  concept  is  important  in  regard  to  our  instru¬ 
ment  as  will  be  discussed  in  connection  with  our  efforts  to  balance  the 
device.  In  the  other  extreme,  for  a  perfectly  balanced  device,  the  model 
assumes  that  the  gradiometer  responds  under  rotations  in  accordance  with 
the  true  physical  gradients.  In  the  intermediate  regime  in  which  the  instru¬ 
ment  is  not  perfectly  balanced  a  systematic  set  of  rotations  about  the 
three  independent  Euler  axes  should  provide,  given  the  magnitude  of  the 
earth's  field,  sufficient  information  to  determine  the  steady  gradient 
field  and  the  magnitude  of  the  gradiometer  imbalance.  The  mathematical 
details  of  this  model  were  presented  in  sections  III  and  VII  of  Podney 


B.  Gradiometer  Response  in  a  Gradient  Field 

To  describe  the  gradiometer  response  during  a  rotation  in  a  gradient 
field  we  must  first  discuss  our  representation  of  a  steady  gradient  field. 
The  gradients  of  a  magnetic  field  in  free  space  in  an  arbitrary  coordinate 
system  can  be  represented  as  a  3  x  3  matrix,  G,  having  elements  g..(r,t) 

i  J 

given  by 

91j(r,t)  *  •$|Xj*l)(r,t)J  (i.j  =  1,  2,  3)  (5) 

where  the  set  of  basis  vectors,  {x^},  is  the  orthogonal  basis  of  the  co¬ 
ordinate  system  and  the  magnetic  field,  b(r,t),  is  a  function  of  the 

position,  r,  and  the  time,  t.  In  this  notation  a  matrix  element,  g.., 

*  J 

represents  the  gradient  in  a  direction  x^  of  the  component  of  magnetic 
field  parallel  to  x..  However,  from  Maxwell's  equations  we  know  that 

s) 

both  the  divergence  and  the  curl  of  a  magnetic  field  in  free  space  are 
zero  so  that  of  the  nine  g^'s  only  five  are  independent.  Specifically 
this  means  that,  in  any  coordinate  system  and  at  each  point  in  space, 
the  matrix  G  is  symmetric  and  has  a  vanishing  trace.  Furthermore,  if 

A 

the  gradients  in  the  basis  { >  are  given  by  the  matrix  G,  then  the  grad- 

A 

ients  in  any  other  arbitrary  basis  {xjl  can  be  found  from  the  matrix 
product: 


G'  =  RGR  (6) 

A 

where  R  represents  the  rotation  which  will  bring  the  basis  {x.j}  into 

-s.  ^ 

coincidence  with  the  basis  {x^},  and  RR  =  RR  =  I  where  the  tilda  denotes 
the  transpose  of  a  matrix  and  1  is  the  identity  matrix.  Now  to  compute 
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the  gradiometer  response  during  a  given  rotation  in  a  known  gradient 
field  we  need  only  define  the  desired  rotation  and  the  coordinate  system 
in  which  the  elements  of  the  gradient  matrix  are  specified. 

First,  let  us  define  a  coordinate  system  fixed  with  respect  to  the 
earth  in  which  the  x  axis  points  west,  the  y  axis  toward  magnetic  north 
and  the  z  axis  downward  and  we  refer  to  the  set  of  three  orthogonal  unit 

A 

vectors  defining  these  directions  as  the  earth  basis,  lx.).  Note  that 

A 

in  this  basis  the  earth's  field,  1^,  has  no  x^  component.  The  second 
coordinate  system  of  interest  is  one  fixed  with  respect  to  the  gradiometer 
with  its  origin  located  midway  between  the  gradiometer  pick-up  loops. 

Here  the  x  axis  lies  in  the  plane  of  the  loops  and  perpendicular  to  the 
line  joining  their  centers,  the  y  axis  lies  perpendicular  to  the  plane 
of  the  loops  and  the  z  axis  is  parallel  to  the  line  joining  the  loops' 
centers  and  points  down  when  the  gradiometer  is  oriented  vertically. 

This  coordinate  system,  specified  by  the  set  of  basis  vectors  {xll,  is 
shown  in  Figure  3  at  some  arbitrary  orientation  with  respect  to  the  earth's 
basis,  {x.l.  Note  that  when  the  gradiometer  is  vertical  with  the  loops 
facing  north  the  above  definition  allows  the  gradiometer  basis  to  coin¬ 
cide  with  the  earth's  basis. 

We  now  ask  what  the  gradiometer  response  would  be  during  a  rotation. 

R,  in  a  steady  gradient  field  if  at  the  start  of  the  rotation  the  gradiom¬ 
eter  basis  were  coincident  with  the  earth  basis.  First  define 

(cos 
sin  <ji 
0 


sin  ij1  0 
cos  ij'  0 
0  1 

(7) 
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EARTH  BASIS 


which  represents  a  rotation  through  an  angle  <p  about  the  x^  axis  of  the 
gradiometer  basis  and 

n  o  o  \ 

X(ip)  =  I  0  cos  ip  -sin  (8) 

\0  sin  cos  ij>  / 

i  » 

which  represents  a  rotation  through  an  angle  \Jj  about  the  x^  axis  of  the 
gradiometer  basis.  We  can  now  represent  the  rotation  which  brings  the 
gradiometer  basis,  (xi),  into  coincidence  with  the  earth's  basis,  {x.} 
by  the  matrix  R(<f>,0,ft)  where 

R(.M,n)  =  z(<j>+tt/2)  x(o)  z(n)  (9) 

u 

where  the  angles  <|>,0  and  ft  are  the  Euler  angles  which  completely  specify 

the  rotation,  R.  Also,  if  the  gradients  in  the  earth's  basis  are  given 

by  the  matrix  G,  the  gradients  in  the  gradiometer  basis,  G',  are  given  J 

by 

G'  =  R(<M,ft)  G  R(4>,0,n)  .  (10) 

But  recall  that  in  the  gradiometer  basis,  the  gradiometer  loops  lie  in 

A  A 

the  (xj.xp  plane  and  the  line  joining  the  centers  of  the  loops  is  paral- 

a  /v 

lei  to  x^.  Hence  in  the  {x.j}  basis,  our  gradiometer  detects  only  the 
gradient  in  the  x^  direction  of  the  component  of  b  along  x£,  that  is  the 
g^  (or,  of  course,  the  g^)  matrix  element.  Consequently  for  a  rotation 
in  a  steady  gradient  field  the  gradiometer  response  as  a  function  of 
<p,0,  and  ft  will  be  precisely  that  found  for  the  g£3  element  of  the  matrix 
G' .  Explicitly  we  find 

* 


(11a) 


r  (0,<!>,n)  =  g^3  =  U( 0 ,<f> )  cos  n  +  V(0 ,4.)  sin  $2 


where 

1  3  1 

u(0 ,4>)  =  (g23cos  <J>-913  sin  4>)cos  20  +  933  +  2^9H"922)C0S  2$  +  912sln  24>] sin  2o 

V(0 ,4>)  =  [g12cos  24»  -  ^(g11-g22)sin  2<|>]sin  0  -  (g13cos  <J>  +  g23sin  4>)cos  0 

(lib) 

Equations  (11a)  and  (lib)  then  represent  our  model  for  the  response  of  a 
perfectly  balanced  gradiometer  during  an  arbitrary  rotation  in  a  steady 
gradient  field  where  the  gradient  elements  g..  are  specified  with  respect 

*  J 

to  the  earth's  basis. 

C .  Gradiometer  Response  From  Imbalance 

From  equation  (4)  we  have  that  the  imbalance  term,  to  first  order, 
has  the  form  $*t)(r,t)  where  the  vector  £  represents  the  assumed  vector 
nature  of  the  gradiometer  imbalance.  As  for  the  gradient  response,  we 
ask  for  the  response  of  the  gradiometer  due  to  imbalance  during  a  rota¬ 
tion,  R,  assuming  that  S  does  not  change  during  the  rotation.  Also  as 
before,  we  will  start  the  rotation  from  the  earth  basis  with  R,  speci¬ 
fied  by  the  Euler  angles  4> ,0 ,n,  as  given  by  equation  (9).  Now  if  R  is 

A 

again  the  rotation  which  takes  the  gradiometer  basis  {xt}  into  the  earth 

A 

basis  {x  .j},  we  have 

t'  =  t  R(<|>,0,i2)  (12) 

and  the  gradiometer  response  from  imbalance,  under  the  rotation  R(<{>,0,i2) 
is  given  by: 
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.1 

rb(o,<j>,n)  =  £’ -6  =  (3r)-S  . 

Explicitly  in  terms  of  4> ,0 ,  and  9  this  becomes, 

rb(0,<J>,fi)  =  0(0,4.)  +  S(0,<f>)cos  9  +  T(e,<j>)sin  9 

where 

0(9 ,<J>)  =  S^tb^cos  0  +  (bjSin  4>  -  b2cos  4>)sin  0] 

S(9 ,<J>)  =  6-|(b-jCos  <J>  +  b2sin  <j>)  +  62[b3sin  0  -  (b^sin  4>  -  b2cos  4>)cos  0] 

T(0,4>)  =  6-|(b3sin  0  -  (b-jSin  <J>  -  b2cos  <t>)cos  0]  -  62(b-|Cos  4»  +  b2sin  <J>) 

(14b) 

where  6^  and  b.  represent  the  individual  components  of  the  vectors  t  and 
S  respectively.  The  components  of  S  in  equations  (14a)  and  (14b)  are 
specified  with  respect  to  the  earth's  basis  while  the  vector  $  is  con¬ 
sidered  fixed  to  the  gradiometer  basis  so  that  the  5^  refer  to  the  gradiom- 
eter  basis,  {x!}.  Note  that  the  expression  for  the  imbalance  response 
contains  terms  in  neither  20  nor  2<f>  as  does  the  expression  for  the  grad¬ 
ient  response 

D.  Total  Gradiometer  Response 

The  total  gradiometer  response  is  now  obtained  trivially  since,  for 
a  linear  device  to  first  order,  it  is  just  the  sum  of  the  gradient  re¬ 
sponse  and  the  imbalance  response.  Summing  equations  (11a)  and  (14a) 
we  find,  the  total  gradiometer  response  to  be  given  by: 


(13) 


(14a) 


r(e,<f>,n)  =  c(o,<j>)  +  A(o,(j))cos  u  +  B(a,(/>)sin  n 
C (0 ,4>)  =  CQ  +  ^[b^cos  0  +  (b^sin  <J>  -  b^cos  )si n  0] 

A(0,<t>)  =  (g23COS  <})-g13sin  4>)cos  20  +  2-{|  923+^(91 1-g22)c°s  2<p  +  g12sin  2i|>]sin  20 
+  6j(b^cos  <f>+b2sin  <)>)  +  62[b^sin  0  +  (b2cos  <{>  -  b^sin  <{»)cos  6] 

B(0,$)  =  [g12cos  2<f>  -  ^-(g11-g22)sin  2$]sin  e  -  (g13cos  <t>  +  g23sin  4>)cos  0 

♦  Cb3s1n  0  -  (bjSin  <*>  -  b2cos  (f>)cos  0]  -  fi^bjCos  <p  +  b2sin  <f>) 

(15) 

Although  we  state  this  result  here  only  for  completeness  we  will  refer  back 
to  this  expression  in  later  sections.  Note  that  with  the  exception  of  the  com¬ 
ponents  of  t ,  the  parameters  on  the  right  hand  side  of  this  equation  are  all 
specified  with  respect  to  the  earth  basis. 


IV.  THE  EQUIVALENT  DIPOLE  MODEL 


In  this  section  we  introduce  the  concept  of  the  "equivalent  dipole" 
for  a  fixed  gradient  field.  As  in  the  previous  section  the  detailed  mathe¬ 
matical  description  of  the  equivalent  dipole  model  is  given  in  Podney  and 
Gillespie,  1977  but  again,  since  the  concept  is  fundamental  to  the  align¬ 
ment  and  noise  suppression  techniques  described  later  we  present 
a  concise  summary  of  the  basic  ideas. 

A.  The  Principal  Axes  of  a  Gradient  Field 

In  addition  to  the  two  coordinate  systems  previously  defined  in  Sec¬ 
tion  III  there  is  a  third  coordinate  system,  which  is  essential  to  the 
concept  of  an  equivalent  dipole  and  is  defined  to  be  coincident  with  the 
principal  axes  of  the  gradient  field  at  a  given  point.  Consequently,  this 
coordinate  system  is  not  defined  in  the  total  absence  of  a  gradient  field. 
However,  since  even  the  unperturbed  magnetic  field  of  the  earth  has  a 
steady  gradient  of  about  10  picotes la/meter  (pT/m)  due  to  its  dipole 
nature  (which  is  detectable  with  our  gradiometer)  we  assume  that,  albeit 
small,  there  always  exists  some  small  gradient  field  with  reference  to  which 
we  can  define  a  set  of  princ<:  1  axes. 

To  define  the  principal  axes  of  a  gradient  field  we  return  to  our 
representation  of  a  gradient  field  in  any  basis  {x..}  as  a  matrix  G  which 
is  real  and  symmetric  with  vanishing  trace.  Also  recall  that  the  grad¬ 
ients  in  any  other  basis  (y^l  can  be  found  from  equation  (6); 

Gy  =  RGR  (6) 
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where  R  is  the  rotation  that  takes  the  {y^}  basis  into  coincidence  with 

a 

the  {Xj}  basis,  and  Gy  has  the  same  properties  as  those  given  above  for 
G.  More  specifically,  for  a  matrix  G  which  is  real  and  symmetric  with 

A 

zero  trace  in  a  basis  {x^}  we  can  always  find  a  real  matrix,  Rp,  such 
that  the  matrix  A  defined  by 


A  =  RpGRp  (16) 

is  diagonal  having  nonzero  diagonal  elements  A^.A^,  and  A^  such  that 
A^+A^+Aj  =  0.  Furthermore,  the  matrix,  Rp,  represents  a  rotation  which 
takes  the  orthogonal  basis  {e..},  in  which  A  is  diagonal,  into  coincidence 
with  the  basis  {x.}.  The  axes  of  the  coordinate  system  in  which  A  is 
diagonal  and  which  are  defined  by  the  set  of  basis  vectors  {e^},  are  re- 
ferred  to  as  principal  axes  of  the  gradient  field.  It  is  important  to  note 
that  as  the  gradients  vary  in  space  from  point  to  point  the  principal  axes 
vary  as  wel 1 . 

Although  conceptually  straightforward,  the  actual  computation  of 
the  principal  axes  and  their  physical  relationship  to  the  initial  basis, 

(x^>  presents  a  complex  geometrical  picture  even  for  a  simple  magnetic 
dipole  field  in  an  ideal  geometry.  However,  since  this  particular  case 
has  special  significance  for  the  present  application  we  will  quote  this  result 
from  Podney  and  Gillespie,  1977.  In  Figure  4  the  dipole  moment  in  is  taken  to 

a 

lie  along  the  x^  axis  and  we  consider  the  gradients  generated  at  a  point  de- 
-►  *  . 

fined  by  a  vector  r  which  lies  in  the  (x^.x^)  plane  and  makes  an  angle  \ 

A  A 

with  the  Xj  axis.  The  gradients  in  the  {x . }  basis  are  then  given  by 
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V  9 


I  cos  x  (1-5  sin2x) 
0 

2 

sin  x(l-5  cos  x) 


0 

cos  x 
0 


2  \ 

sin  x(l-5  cos  x) 

0 

cos  x(‘2+5  sin2x)j 


(17) 


where 


Mo  /3m, 

9  '  • 
r 


(18) 


Here  uQ  is  the  permeability  of  free  space,  r  =  |r|,  and  m  =  |m|.  The 
eigenvalues  of  this  matrix  are  found  by  obtaining  the  roots  of  the  equation 


det(Gd-AI)  =  0  (19) 

where  A  is  a  parameter  and  I  is  the  identity  matrix.  As  shown  previously 
(Podney  and  Gillespie,  1977)  the  roots,  A^.A^.  and  A^  are: 

Aj  =  (g/2) 1(4+5  cos2x)2  -  cos  x) 

A2  =  g  cos  x 

A3  =  -(g/2) [(4+5  cos2x)^  +  cos  x)  (20) 


Hence  the  diagonal  matrix  A  having  elements  A-j ,  A2  and  A^  gives  the  elements 
of  the  gradient  matrix  G  in  the  principal  axes  basis,  {e^}.  For  this 

A 

choice  of  the  basis  {x^}.  Figure  4  shows  the  principal  axes  for  a  con¬ 
stant  value  of  | r |  at  several  values  of  x  where  0  <  x  £ir.  The  details 
of  this  derivation  can  be  found  in  Appendix  A  of  Podney  and  Gillespie, 

1977. 
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B.  The  Equivalent  Dipole  of  a  Gradient  Field 

Simply  stated  the  concept  of  the  "equivalent  dipole  of  a  gradient 
field"  means  that  at  any  point  in  free  space  (where  div  S  =  curl  S  =  0) 
the  gradient  field  is  identical  to  the  gradients  produced 
by  some  dipole  if  the  dipole  position,  orientation  and  strength  are  properly 
chosen.  To  illustrate  this  consider  an  arbitrary  gradient  field  at  a  point 
in  free  space  which  has  a  gradient  field  specified  by  the  five  independent 
elements  of  the  matrix  G'.  Let  us  now  determine  the  eigenvalues  of  G',  .1 

Aj ,  and  order  them  so  that  Aj  >  A£  —  The  9radient  field  can 

now  be  specii .ed  in  the  principal  axes  by  the  diagonal  matrix,  A',  having 

nonzero  diagonal  elements  A^,A£,A^  and  we  would  like  to  show  that  for  the  approp- 
riate  choice  of  the  parameter, x  .  the  matrix  A'  is  exactly  equivalent 
to  the  diagonal  matrix  A  which  represents  the  gradient  field  from  a  dipole. 

In  the  first  place  we  know  that  t 


X-j  +  *2  +  *3  =  *1  i  *2  —  *3 

(21a) 

A-j  +  X2  +  X^  =  0,  Aj  >  >  A3 

(21b) 

* 

so  that  we  must  have 

A]»  Aj  >  0  Xj>  Ag  <  0 

(22) 

• 

and  X2  and  A£  can  be  either  less  than  or  greater  than  zero.  Let  us  first 
assume  A£  >  0  and  consider  the  ratios  X-j/X^  and  A-j/A^.  From  equations 
(21)  and  (22)  we  have 

X-j/X2  >  1  *  ^3/^2  —  >  0  (23) 
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(24) 


and  we  now  look  for  the  ratio  X-j/A^  such  that 

Aj/A£  =  A/A2  . 

From  the  expressions  for  A-j ,  A^.  and  A^  in  equations  (20)  we  have  plotted 
In  Figure  5  Aj/A2  (solid  lines)  and  A3/A2  (dashed  lines)  as  functions 
of  x-  By  the  same  arguments  leading  to  equation  (23),  for  A^  >  0  we 
have 


A1/A2  >  1  ,  A3/A2  <  -2  ;  A2  >  0  (25) 

Since  A-j/A2  and  A|/A2  vary  over  exactly  the  same  interval  for  A2»  A£  >  0, 
for  any  possible  Aj  and  A£  we  can  always  find  a  value  of  A^/A2  (and  x) 
such  that  equation  (24)  is  satisfied;  which  also  requires  from  equations 
(21a)  and  (21b)  that 


VX2  =  X3/X2  (26) 

A  trivial  manipulation  of  equations  (21a),  (21b),  (24),  and  (26)  then 
shows  that 


A 

A 


1 

I 

1 


g' 


(27) 


or 


A  =  g'A'  (28) 

where  the  factor  g'  just  represents  a  scaling  factor  given  by  the  magnitude 
of  the  magnetic  dipole  moment.  An  exactly  analogous  argument  can  be  applied 
to  the  case  where  A£  <  0,  and  if  A£  =  0  we  have  Aj  =  A3  and  similarly 
choose  A2  =  0  so  that  A-j  =  -A3  to  again  arrive  at  equation  (28).  Hence 
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Figure  5.  Graph  of  Xi/X2  and  X3/X2  as  a  function  of  x  where  Xi,X2,  and  Xj  are  the  eigenvalues  of  the  grad¬ 
ient  matrix,  Gj  where  G^  represents  the  gradients  in  a  basis  {x.}  from  a  magnetic  dipole,  ifi,  where 
m  lies  along  x3.  The  angle  x  is  the  angle  between  the  position  vector  t  and  the  dipole  moment  ifi. 


as  a  consequence  only  of  Maxwell's  equations  for  free  space  which  requires 
div  b  =  curl  b  =  0,  and  the  ordering  of  the  eigenvalues  which  constitute  A 
and  A',  we  see  that  the  gradient  field  at  any  point  in  free  space  is  exactly 
equivalent  to  a  dipole  field  if  the  proper  values  of  x  and  g'  are  chosen. 

In  practice,  when  the  calculation  is  performed  we  find  that  the  orientation 
of  the  orthogonal  set  of  principal  axes  is  unique.  However,  because  the  def¬ 
inition  of  a  right  handed  coordinate  system  with  this  orientation  is  arbit¬ 
rary  with  respect  to  90°  rotations  we  find  that  there  are  actually  four  pos¬ 
sible  locations  and  orientations  for  the  equivalent  dipole  of  the  gradients 
at  any  particular  point.  Nonetheless  in  regard  to  the  gradients  at  the  point 
of  interest,  any  of  the  four  possible  equivalent  dipoles  is  suitable  and  we 
merely  choose  that  which  has  the  most  convenient  geometrical  location. 

This  completes  our  presentation  of  the  basic  theoretical  concepts  used 
in  our  balancing,  alignment  and  noise  suppression  experiments  and  in  the  re¬ 
maining  sections  we  will  refer  freely  to  these  ideas. 


■  I 
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V.  NOISE  SOURCES  AND  NOISE  SUPPRESSION  TECHNIQUES 


During  the  conceptual  phases  of  the  experiment  we  had  recognized  that 
a  major  difficulty  in  performing  any  magnetic  measurements  from  the  NOSC 
oceanographic  tower  would  be  the  immediate  presence  of  the  steel  structure 
of  the  tower  itself  and  the  resultant  magnetic  moment.  Since  ionospheric 
disturbances  of  the  earth's  magnetic  field  occur  at  the  frequencies  expec¬ 
ted  for  internal  waves  the  fluctuations  of  the  magnetic  moment  of  the  tower, 
driven  by  the  earth's  field,  eliminated  the  possibility  of  performing 
measurements  directly  from  the  tower  proper.  However,  preliminary  measure¬ 
ments  of  the  magnetic  and  gradient  fields  of  the  tower  indicated  that 
at  a  distance  of  about  ?5  meters  from  the  tower's  center,  the  fields  and 
gradients  of  the  tower  conformed  reasonably  well  to  a  simple  model  consist¬ 
ing  of  a  dipole  in  an  otherwise  uniform  magnetic  field.  From  these  measure¬ 
ments  we  then  deduced  the  location,  orientation  and  magnitude  of  the  equiv¬ 
alent  dipole  of  the  tower  and  on  the  basis  of  this  information  we  developed 
the  noise  suppression  techniques  described  earlier  (Gillespie  and  Podney, 
1976).  The  actual  experiment  was  performed  from  a  cantilevered  structure 
attached  to  the  south  face  of  the  NOSC  tower  which  supported  the  gradio- 
meter  on  a  platform  approximately  25  meters  from  the  center  of  the  tower 
and  7  meters  above  the  water.  The  preliminary  magnetic  measurements  at 
the  tower  and  the  cantilevered  boom  have  been  described  in  more  detail 
previously  (Podney  and  Gillespie,  1977).  In  this  section  we  will  delineate 
the  various  noise  sources  we  anticipated  with  this  geometrical  arrangement 
and  the  techniques  used  to  suppress  the  noise  from  each  source. 

A.  Noise  Sources 

The  major  sources  of  noise  at  frequencies  below  5  milliHertz  (mHz)  with 
experimental  configuration  can  be  categorized  into  four  basic  sources: 


(1)  Magnetization  currents  by  which  we  mean  the  fluctuations  of  the 
permanent  dipole  moment  of  the  steel  tower  in  response  to  changes  in  the 
spatially  uniform  field  of  the  earth. 

(2)  Conduction  currents  or  eddy  currents  which  flow  through  any 
electrically  conducting  material  and  which  are  also  driven  by  fluctuations 
of  the  earth's  field. 

(3)  Relative  motion  between  the  gradiometer  and  the  magnetic  dipole 
moment  of  the  tower  as  a  result  of  small  oscillations  or  rotations  of  the 
cantilevered  boom.  In  our  configuration  this  represented  the  most  severe 
source  of  noise  due  to  the  large  ambient  gradient  of  the  tower  (approxi¬ 
mately  200  nT/m). 

(4)  Gradiometer  imbalance  which  produces  a  signal  at  the  gradiometer 
output  during  small  rotations  of  the  gradiometer  in  the  uniform  field 

of  the  earth.  The  imbalance  of  the  gradiometer  arises  from  small  differ¬ 
ences  in  the  area  and  orientation  of  the  paired  pick-up  loops  resulting 
in  a  magnetometer- like  response  during  rotations  in  a  spatially  uniform 
field.  Thus  noise  can  arise  in  an  unbalanced  gradiometer  from  either  ro¬ 
tations  of  the  gradiometer  in  a  spatially  uniform  field  or  from  time  depen¬ 
dent  fluctuations  of  a  spatially  uniform  field. 

The  anticipated  internal  wave  spectrum  lies  at  most  a  factor  of  about 
10  above  the  instrument  noise  at  frequencies  of  1  to  5  mHz  (Rodney  and  Gil¬ 
lespie,  1977).  Hence,  to  assure  a  reasonable  probability  of  detecting  the 
fluctuating  magnetic  gradients  from  internal  waves  it  was  necessary  to  reduce 
the  noise  contributions  from  the  various  sources  to  at  least  a  level  compar¬ 
able  to  that  of  the  SQUID  sensor  itself. 


B.  Noise  Suppression  Techniques 

To  suppress  the  noise  from  the  external  sources  outlined  above  we 
used  the  concept  of  the  equivalent  dipole  introduced  in  the  last  section 
in  which  we  represent  the  gradients  of  a  magnetic  field  at  a  point  in  free 
space  as  equivalent  to  those  produced  by  an  appropriate  magnetic  dipole. 

Thus  for  our  experimental  configuration  in  which  the  gradiometer  is  25 
meters  from  the  center  of  the  tower  we  modeled  the  ambient  magnetic  field 
and  its  gradients  as  a  spatially  uniform  field  plus  a  dipole  contribu¬ 
tion  from  the  tower  itself.  To  first  order  this  provided  a  workable  pic¬ 
ture  with  which  to  approach  the  problem  of  noise  suppression.  To  mini¬ 
mize  the  effects  of  eddy  currents,  particularly  very  near  the  gradiometer, 
and  to  maintain  a  magnetically  clean  environment  the  outermost  13  meters 
of  the  cantilevered  boom  were  constructed  entirely  of  nonmagnetic,  non¬ 
conducting  materials  consisting  primarily  of  fiberglass  and  fiberglass 
resins.  The  only  eddy  current  contributions  then  were  those  produced  in 
the  tower  itself  and  in  the  aluminum  portion  of  the  cantilevered  boom  by 
fluctuations  in  the  earth's  field.  Since  the  typical  fluctuations  in  the 
earth's  field  (except  during  periods  of  high  solar  activity)  are  only  of 
order  10  nT,  or  about  .02  percent,  we  represent  the  magnetic  gradients 
from  magnetization  currents  and  eddy  currents  in  the  tower  as  small  changes 
in  the  position,  orientation,  and  magnitude  of  the  steady  tower  dipole.  If  we 
also  include  in  this  picture  small  motions  of  the  gradiometer  due  to  oscil¬ 
lations  or  twisting  of  the  cantilevered  boom  then  all  of  the  above  noise 
contributions  except  those  from  imbalance  can  be  considered  as  small  per¬ 
turbations  on  the  steady  gradient  field  from  the  tower.  Furthermore, 
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such  perturbations  can  be  represented  as  small  fluctuations  in  the  posi¬ 
tion,  orientation,  and  strength  of  the  equivalent  dipole. 

Justification  for  this  representation  is  provided  by  the  following 
considerations.  Clearly  the  permanent  magnetic  moment  of  the  tower  which 
arises  from  the  stationary  steel  structure  residing  in  the  earth's  field 
can  be  expected  to  vary  only  on  the  same  scale  as  does  the  earth's  field 
itself;  that  is,  of  order  .02  percent.  Moreover,  the  stability  of  our 
cantilevered  boom  is  better  than  10~4  radians  (RMS)  in  the  frequency  range 
of  interest  which  will  produce  an  apparent  gradient  fluctuation  of  order 
.01  percent.  Hence  these  effects  can  certainly  be  represented  as  small 
fluctuations  in  the  tower  gradients. 

Gradient  fluctuations  arising  from  eddy  currents  however  will  depend 
on  various  geometrical  factors.  We  first  consider  an  eddy  current  path 
in  the  aluminum  portion  of  the  cantilevered  boom  closest  to  the  gradiom- 
eter.  The  voltage  induced  in  a  current  loop  will  be 

V  =  ~  =  k6A  (29) 

where  B  is  the  fluctuating  field,  w  is  the  angular  frequency  of  the  fluc¬ 
tuations  in  radians/sec  and  A  is  the  area  of  the  current  loop.  The  cur¬ 
rent  induced  is  then 


V  _  u8A 
R  "  R 


(29a) 


where  R  is  the  total  resistance  around  the  loop.  The  magnetic  gradients 
at  a  distance,  r,  from  the  loop  can  now  be  approximated  by  considering 
the  current  loop  to  be  a  dipole  source: 
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where  the  dipole  moment  m  is  just 


m  =  iA  = 


To  describe  some  average  current  path  in  the  aluminum  portion  of  the  boom 


we  selected  the  following  values: 


radius  of  loop  =  5  meters 


B  =  10  nT 


w  =  10"2  Hz 


r  =  15  meters 


R  =  10"3n 


The  value  for  the  resistance,  R,  was  arrived  at  by  considering  the  tubu¬ 


lar  aluminum  in  the  boom  to  have  a  V  wall  and  a  resistivity  of  5x10  fi-cm*. 


These  values  give  a  total  gradient  of  about  3x1 0"3  nT/m  at  the  gradiom- 


eter.  Gradients  from  eddy  currents  in  the  tower  would  be  even  smaller 


since,  although  the  current  loop  area  may  be  larger,  the  distance  from 


the  gradiometer  is  greater  and  the  resistance  of  the  steel  is  about  a 


factor  of  10  greater  than  that  of  the  aluminum.  Since  these  effects  are 


of  roughly  the  same  order  of  magnitude  as  the  fluctuations  due  to  magnet¬ 


ization  currents  and  boom  stability,  we  treated  them  all  as  small  pertur¬ 


bations  on  the  steady  gradient  field  of  the  tower. 
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Two  approaches  can  now  be  used  to  reduce  the  noise  generated  in  the 
gradiometer  by  these  effects.  First,  by  determining  the  position,  orien¬ 
tation  and  strength  of  the  equivalent  dipole  it  might  be  possible  to 
position  a  set  of  nested  coils  with  which  to  cancel,  at  the  site  of  the 
gradiometer,  the  steady  gradients  from  the  tower.  This  would  produce  a 
region  near  the  gradiometer  with  a  reduced  gradient  field.  Since  the 
gradiometer  response  to  small  angle  rotations  in  a  gradient  field  is  pro¬ 
portional  to  the  product  of  the  gradient  and  the  angular  displacement, 
the  noise  produced  in  the  gradiometer  as  a  result  of  motions  of  the  boom 
would  also  be  reduced. 

In  addition  to  this  technique  it  is  also  possible  to  specify  certain 
orientations  of  the  gradiometer  with  respect  to  the  equivalent  dipole  such 
that  the  gradiometer  response  to  small  perturbations  of  the  dipole  posi¬ 
tion  and  orientation  is  minimized.  To  picture  this  let  us  consider  the  gradiom 

A 

eter  basis  defined  in  Section  III  (such  that  is  parallel  to  the  line  join- 

A 

ing  the  pick-up  loops  and  xj  is  in  the  plane  of  the  pick-up  loops).  Then, 
on  the  basis  of  previous  calculations  (Podney  and  Gillespie,  1977),  we 
expected  this  so-called  Minimum  Noise  Orientation  (MNO)  to  be  defined  as 
that  orientation  of  the  gradiometer  such  that  the  vector  from  the  gradiom- 

/v 

eter  to  the  position  of  the  dipole  be  parallel  to  x|  and  the  gradiometer  pick- 
up  loops  lie  in  the  plane  of  x^'  and  the  dipole  moment,  m,  of  the  equiva¬ 
lent  dipole.  Furthermore,  in  this  orientation  the  gradiometer  is  tota'My 
insensitive  to  changes  in  the  dipole  orientation  so  that  small  changes 
in  the  orientation  of  the  equivalent  dipole  of  the  tower  produced  by 


magnetization  currents  or  eddy  currents  should  produce  no  response  at  all  in 
the  gradiometer. 

As  we  discuss  in  detail  in  Section  VI,  the  gradiometer  orientation  correspond¬ 
ing  to  the  M.N.O.  proved  to  be  somewhat  different  than  originally  anticipated. 
Nonetheless,  the  concept  proved  to  be  a  powerful  technique  with  which  to  reduce 
the  noise  generated  by  the  presence  of  the  tower  and  in  the  final  phases  of  the 
experiment  we  depended  on  it  exclusively. 

To  conclude  this  discussion  of  noise  suppression  techniques  we  must  ad¬ 
dress  the  remaining  question  of  noise  introduced  into  the  gradiometer  by  im¬ 
balance.  For  typical  fluctuations  of  the  earth's  field  of  about  10  nT  at  the 
frequencies  of  interest  an  imbalance  of  4x10"^  m~  ^  (or  1  part  in  10^)  produces 
an  effective  signal  from  the  gradiometer  of  about  4xl0"4  nT/m  or  about  .4  pT/m. 

This  is  at  least  a  factor  of  5  below  the  SQUID  sensor  noise  at  frequencies  of 
1-5  mHz.  Only  a  modest  level  of  balance  then  is  necessary  to  eliminate  noise 
from  fluctuations  of  the  earth's  field.  A  more  serious  problem  concerning 
the  balance  is  the  effect  of  small  rotations  of  the  gradiometer  in  the  earth's 
field  of  4.5xl04  nT.  In  this  case,  the  gradiometer  signal  will  be  given  approx¬ 
imately  by 

ig  =  b  6  A0  (32) 

o 

where  bQ  is  the  earth's  field,  6  is  the  magnitude  of  the  imbalance  vector 

and  A0  represents  a  small  rotation  of  the  gradiometer  with  respect  to  bQ. 
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Again  for  an  imbalance  of  6  =  4x10  m  and  bQ  =  4.5x10  nT  we  have 

Ag  =  (2  nT/m)A0  .  (33) 

The  stability  of  our  cantilevered  boom  at  the  frequencies  of  interest  is  about 


-5  ? 

3x10  radians  giving  an  approximate  equivalent  noise  of  Ag  <  6x10  pT/m, 

again  well  below  the  instrument  noise  at  an  RMS  level  of  about  2  pT/m.  Con- 

-5  -1 

sequently,  for  our  experiment  a  rather  modest  balance  of  4x10  m  was  suf¬ 
ficient  to  insure  that  the  noise  in  the  gradiometer  from  imbalance  effects 
was  well  below  the  instrument  noise  level.  A  more  complete  description  of 
our  efforts  to  balance  the  instrument  and  of  some  unexplained  effects  we  ob¬ 
served  during  balancing  of  the  gradiometer  will  be  presented  in  a  later  sec¬ 
tion. 
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VI.  EXPERIMENTAL  NOISE  SUPPRESSION  RESULTS 


As  outlined  in  the  previous  section,  the  basic  approaches  to  the  noise 
suppression  problem  fall  into  three  general  categories:  (1)  gradiometer  bal¬ 
ance,  (2)  alignment  of  the  gradiometer  with  respect  to  the  steady  gradient 
field,  and  (3)  cancellation  of  the  steady  gradients  using  a  set  of  three  nested 
coils.  For  our  instrument  the  behavior  of  the  balance  proved  to  be  a  rather 
complicated  phenomenon  which  we  describe  in  detail  in  Section  VII.  Nonethe¬ 
less,  we  achieved  and  maintained  a  balance  of  order  4x10"^  m"^  which  was  suf¬ 
ficient  to  insure  that  the  noise  introduced  into  the  gradiometer  from  imbalance 
was  well  below  the  inherent  sensor  noise  at  the  frequencies  of  interest. 

In  light  of  the  above,  the  discussion  in  this  section  will  deal  only 
with  gradiometer  alignment  and  cancellation  of  the  steady  gradients  near  the 
gradiometer.  The  precision  with  which  gradient  cancellation  can  be  implemented 
is  directly  determined  by  the  precision  with  which  we  can  measure  the  steady 
gradients  and  determine  the  parameters  specifying  their  equivalent  dipole. 

Our  first  approach  to  this  problem  consisted  of  performing  a  systematic  set  of 
rotations  with  the  gradiometer,  analyzing  ithe  gradiometer  output  in  the  con¬ 
text  of  equation  (15)  which  represents  our  model  for  the  gradiometer  response, 
and  computing  the  approximate  values  of  the  ambient  gradients.  Then  starting 
from  this  first  estimate  a  more  precise  measurement  of  the  gradients  would  be 
performed  using  an  iterative  technique,  again  predicated  upon  the  use  of  various 
gradiometer  rotations.  In  short,  the  initial  set  of  gradiometer  rotations 
did  indeed  provide  a  reasonable  first  estimate  of  the  gradient  field  but  the 
iterative  technique  totally  failed  to  provide  the  precision  required  to  compute 
the  gradiometer  and  coil  alignments  with  sufficient  accuracy.  As  a  result 
of  this  failure  we  subsequently  developed  a  somewhat  more  empirical  technique 
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to  determine  the  optimum  alignment  of  the  gradiometer  in  a  gradient  field. 

Tests  at  La  Posta  demonstrated  the  effectiveness  of  the  technique  and  dur¬ 
ing  the  tower  experiments  we  found  that  proper  alignment  alone  was  sufficient 
to  suppress  all  of  the  noise  from  external  sources  to  a  level  comparable  to 
that  of  the  inherent  gradiometer  noise.  In  this  section  we  will  describe 
this  noise  suppression  technique  and  present  the  results  of  some  of  the  empir¬ 
ical  measurements  used  to  determine  the  so-called  Minimum  Noise  Orientation 
(MNO)  and  some  comparative  gradiometer  noise  spectra  demonstrating  the  effect¬ 
iveness  of  noise  suppression. 

A.  Concept  o f  the  Minimum  Noise  Orienta ti on 

The  basic  concept  of  the  MNO  is  that  there  exists  some  orientation  of 
the  gradiometer  with  respect  to  the  ambient  gradient  field  at  which  the  grad¬ 
iometer  response  to  an  arbitrary  small  angle  rotation  will  be  minimized. 
Mathematically  the  MNO  can  be  found  by  calculating  the  derivative  of  the  grad¬ 
iometer  response  with  respect  to  each  of  the  five  parameters  specifying  the 
equivalent  dipole  of  the  gradient  field  and  finding  the  three  Euler  angles 
at  which  all  five  derivatives  are  zero.  In  fact,  there  is  no  position  at 
which  all  five  of  the  derivatives  vanish  simultaneously  (as  discussed  in 
Podney  and  Gillespie,  1977)  although  there  does  exist  a  position  at  which  four 
of  the  derivatives  can  be  made  to  disappear.  This  is  precisely  the  orienta¬ 
tion  described  in  Section  V  in  which  the  position  of  the  equivalent  dipole 

A 

lies  on  the  x^  axis  of  the  gradiometer  basis  and  the  gradiometer  pick-up  loops 

*  -> 

lie  in  the  plane  defined  by  xj  and  the  dipole  moment,  m,  of  the  equivalent 
dipole.  At  this  orientation  the  gradiometer  is  totally  insensitive  to  changes 
in  the  dipole  orientation  and  strength  (regardless  of  their  magnitude)  and 


is  also  insensitive  to  small  changes  in  the  dipole's  position  within  the  plane 


of  the  loops.  Motions  of  the  dipole  out  of  this  plane  however  produce  a  grad- 
ioineter  response  proportional  to  the  angular  displacement  of  the  dipole  from 
the  plane  with  respect  to  the  position  of  the  gradioineter. 

It  should  be  noted  here  that  since  all  five  of  the  derivatives  do  not 
vanish,  attention  must  be  given  to  the  relative  magnitudes  of  the  derivatives 
as  well  as  the  predominant  modes  of  relative  motion  between  the  gradiometer 
and  the  equivalent  dipole.  For  example,  in  a  simple  picture  if  the  relative 
positions  of  the  gradiometer  and  equivalent  dipole  are  rigidly  maintained  but 
the  orientation  and/or  magnitude  of  the  dipole  undergoes  large  fluctuations, 
the  orientation  described  above  is  ideal  due  to  the  gradiometer's  total  in¬ 
sensitivity  to  the  dipole  orientation  and  strength.  On  the  other  hand  if  the 
gradient  field  and  its  equivalent  dipole  are  relatively  constant  but  the  grad¬ 
iometer  moves  significantly  with  respect  to  the  gradient  field,  it  may  be  ad¬ 
vantageous  to  use  an  orientation  where  the  gradiometer  is  insensitive  to  motions 
out  of  the  plane  although  two  other  derivatives  may  be  non-zero  This  point 
has  particular  relevance  to  our  measurements  as  we  will  discuss. 
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B. 


Determination  of  the  Minimum  Noise  Orientation 


As  mentioned  above,  experimental  limitations  prevented  us  from  measure- 
# 

ing  the  steady  gradient  field  at  the  position  of  the  gradiometer  with  suf¬ 
ficient  precision  to  calculate  t.  .  appropriate  MNO.  An  alternative  ap¬ 
proach  was  then  suggested  by  the  very  nature  of  the  problem.  Specifically, 
noise  is  introduced  into  the  gradiometer  as  a  result  of  relative  motions 
between  the  gradiometer  and  the  equivalent  dipole.  Furthermore,  for  small 
relative  motions,  aligning  the  gradiometer  to  the  MNO  will  minimize  the 
noise  produced  by  such  motions.  Since  this  concept  should  depend  only 
on  relative  orientations  and  motions  and  be  essentially  independent  of 
frequency  (as  long  as  the  gradiometer  has  the  same  average  position  and 
orientation  for  each  frequency  component)  it  should  be  possible  to  moni¬ 
tor  the  noise  level  of  the  gradiometer  at  a  few  Hertz  as  a  function  of  the 
gradiometer  orientation. 

Using  a  Spectral  Dynamics  Model  330  real-time  spectrum  analyzer 
to  measure  the  gradiometer  noise  spectrum  from  0.1  to  25  Hz,  our  measure¬ 
ments  were  performed  as  follows.  Starting  at  the  MNO  estimated  from  initial 
rough  measurements  of  the  steady  gradient  field,  a  gradiometer  noise  spectrum 
is  made  at  incremental  steps  of  one  of  the  Euler  angles,  say  $,  while  hoi  'ing  G 
and  ft  constant.  The  noise  level  at  some  frequency  or  over  some  frequency  range 
is  then  measured  for  each  value  of  4>  and  plotted  versus  4>  (since  we  expect  the 
noise  to  first  order  to  be  proportional  to  the  angular  misalignment).  The 
process  is  then  repeated  for  the  other  Euler  angles  0  and  ft,  and  subsequent 
iterations  eliminate  possible  interactions  among  the  three  angles. 
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Figure  6  shovis  results  of  the  final  iteration  for  one  of  our 
first  tests  of  this  technique  at  the  La  Posta  site  using  an  iron  sphere 
five  feet  in  diameter  positioned  approximately  12  feet  away  from  the 
gradiometer.  To  provide  a  more  appropriate  simulation  of  the  NOSC  tower 
environment  the  iron  sphere  was  cushioned  with  foam  rubber  to  provide  a 
motional  aspect  (driven  by  the  wind)  with  respect  to  the  gradiometer. 

The  very  strong  dependence  of  the  noise  level  on  gradiometer  alignment 
is  clear  from  Figure  6.  Using  the  MNO  determined  in  this  manner  we  sub¬ 
sequently  achieved  a  gradiometer  noise  level  in  the  presence  of  the  large 
gradient  field  of  the  iron  sphere  essentially  equivalent  to  the  inherent 
gradiometer  noise  level  in  the  range  of  1  to  10  mHz.  In  Figure  7 
the  solid  line  shows  the  inherent  gradiometer  noise  spectrum  made  in  a 
magnetically  clean  environment  (after  instrument  modifications  were 
completed)  while  the  dashed  line  shows  the  noise  spectrum  with  the  iron 
sphere  at  a  distance  of  12  feet  and  with  the  gradiometer  at  the  MNO.  At 
this  orientation  the  gradiometer  pick-up  loops  were  essentially  edge-on  to¬ 
ward  the  iron  sphere.  For  comparison,  the  noise  level  with  the  pick-up  loops 
facing  the  iron  sphere  is  shown  in  Figure  7  by  the  broken  line. 

C.  Minimum  Noise  Orientation  of  the  NOSC  Tower 

The  data  collected  at  the  La  Posta  site  were  presented  primarily 
to  demonstrate  the  remarkable  effectiveness  of  this  technique  even 
during  our  initial  tests.  Nonetheless  the  La  Posta  experiments 
provided  only  a  proving  ground  and  the  object  of  the  tower  experiment 
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POWER  SPECTRAL  3ENS.7Y  [(pT/nn)2/H2: 


n<jurc  /.  noise  spectra  comparing  no  I  so  levels  with  yrailiometo 

at  Minimum  No  iso  Orientation  (  )  .uni  with  loops  f.u'imi  iron 

sphere  (  -  ).  Solid  lino  is  instrument  noise  in  mag¬ 

netically  clean  environment. 


in  May  of  1978  was  to  demonstrate  the  feasibility  of  i nip lemen tiny  the  same 
methodology  at  the  tower  location  and  to  achieve,  in  the  absence  of  internal 
waves,  a  gradiometer  noise  level  on  the  tower  comparable  to  the  inherent  in¬ 
strument  noise.  In  fact,  the  technique  did  work  well  and  the  data  from  the 
tower  experiment  showing  the  MNO  measurement  and  the  accompanying  comparative 
gradiometer  noise  spectra  are  shown  in  Figures  8  and  9.  Here  again  the  MNO 
measurements  represent  the  final  iteration  and  give  the  values  for  0Q,  $o, 
and  Qq  that  we  used  throughout  the  remainder  of  the  May  measurements  and  all 
of  the  internal  wave  measurements  during  August.  In  Figure  9  the  solid  line 
denotes  the  instrument  noise  in  a  magnetically  clean  environment,  the  dashed 
line  shows  the  noise  with  the  gradiometer  at  the  MNO  (with  the  loops  essentially 
edge  on  to  the  tower),  and  the  broken  line  provides  a  comparative  spectrum  in 
which  the  gradiometer  pick  up  loops  were  facing  the  tower. 

A  final  comment  is  required  on  the  values  of  0O,  <t>Q,  and  determined 
for  the  MNO  by  the  above  measurement.  If  the  equivalent  dipole  of  a  gradient 
field  is  known  we  can  calculate  the  gradiometer  orientation  at  which  we  ex¬ 
pect  the  lowest  gradiometer  sensitivity  to  fluctuations  in  the  dipole.  In 
particular  we  look  for  an  orientation  at  which  the  derivative  of  the  gradiom¬ 
eter  output  with  respect  to  each  of  the  five  independent  dipole  parameters 
vanishes.  Since  there  is  no  orientation  at  which  all  five  of  the  derivatives 
vanish  simultaneously  we  initially  guessed  that  the  MNO  would  occur  at  the 
orientation  where  four  of  the  derivatives  are  zero.  Using  a  technique  which 
we  describe  in  Section  VIII  ,  we  measured  the  ambient  gradient  field  with  our 
gradiometer  and  computed  its  equivalent  dipole  to  an  accuracy  corresponding  to 
an  error  of  about  three  degrees  in  position  and  orientation.  Based  on  this 
computation  we  calculated  an  MNO  for  the  tower  gradients  given  by  0o  =  25.9° 
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MINIMUM  NOISE  ORIENTATION 


6  (DEG)  <P(DEG)  R(DEG) 

Figure  3.  Gradlometer  noise  as  a  function  of  gradiometer  orientation  at  NOSC  tower.  In  each  graph  two  angles 
were  held  constant  while  the  third  angle  was  varied.  Solid  lines  indicate  linear  fit  to  data  and 
arrows  indicate  values  determined  for  minimun  noise  orientation. 
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Figure  9.  Gradiometer  noise  spectra  comparing  noise  levels  with  gradiometer 

at  minimun  noise  orientation  ( - )  and  with  loops  facing  tower 

( - ).  Solid  line  is  instrument  noise  in  a  magnetically 

clean  environment. 


<J>0  =  198°  and  i2Q  =  74.8°,  values  which  are  significantly  different  from  those 
determined  from  the  measurements  displayed  in  Figure  8.  Although  we  estimate 
the  accuracy  of  the  calculated  values  to  be  of  order  +3°  and  the  accuracy  of 
the  measured  values  from  Figure  8  to  be  of  order  ±1°,  the  possible  errors  do 
not  account  for  the  differences,  particularly  the  difference  in  tilt  angle. 

Although  the  MNO  noise  measurement  was  repeated  several  times  (even  starting 
with  values  of  0  near  20°),  we  consistently  found  a  well  defined  minimum  in 
the  noise  spectrum  at  values  of  0  near  8°.  The  discrepancy,  which  lies  in  the 
criterion  we  initially  used  to  guess  the  MNO,  is  resolved  in  Section  IX. 

D.  Gradient  Cancellation  Experiments 

Although  we  found  during  field  tests  at  La  Posta  that  alignment  by  it¬ 
self  could  provide  sufficient  noise  suppression,  cancellation  of  the  gradient 
field  was  still  a  desirable  goal  from  the  standpoint  of  balancing  the  instru¬ 
ment  when  in  the  presence  of  a  large  gradient.  Furthermore,  a  reduced  gradient 
field  should  make  the  gradiometer  less  sensitive  to  small  rotations  in  the 
gradient  field  and  require  less  precision  when  performing  the  gradiometer 
alignment.  During  our  work  at  the  tower  during  May,  1978  we  experimented  with 
the  coil  alignment  and  gradient  cancellation  with  regard  to  balancing  the  instru¬ 
ment  in  the  presence  of  the  tower  gradients.  Measurements  with  the  gradiometer 
indicated  residual  gradients  of  order  12  nT/m  with  the  coils  adjusted  to  can¬ 
cel  the  tower  gradients.  For  this  coil  configuration  we  collected  two  grad¬ 
iometer  spectra  for  comparison;  one  with  the  coils  off  and  another  with  the 
coils  on  keeping  all  other  parameters  the  same.  For  these  data  the  gradiom¬ 
eter  was  aligned  to  the  calculated  MNO  based  on  the  position  and  orientation 
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of  the  coil  dipole  such  that  the  center  of  the  coils  lay  on  the  xj  axis  of  the 
gradiometer  basis  with  the  gradiometer  pick-up  loops  in  the  plane  defined  by 

A 

x^  and  the  coil  dipole  moment,  mc. 

^  i 
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These  data  are  presented  in  Figure  10  in  which  the  coils  off  data  is 
shown  by  the  dashed  line,  the  coils  on  data  by  the  broken  line  and  the  solid 
line  is  the  gradioineter  noise  observed  on  the  tower  at  the  experimentally  de¬ 
termined  MNO  (from  Figure  9).  Comparison  of  the  dashed  and  broken  curves  shows 
the  noise  reduction  achieved  with  the  coils  to  be  only  about  a  factor  of  3 
instead  of  the  expected  factor  of  15.  The  discrepancy  probably  arises  from 
relative  motions  between  coils,  gradiometer,  and  tower  and  from  fluctuations 
in  the  gradient  field  of  the  coils.  Since  the  coils  are  located  on  the  canti¬ 
levered  boom,  twistinq  and  oscillations  of  the  boom  produce  relative  motions 
of  the  coils  with  respect  to  both  tower  and  gradiometer.  In  this  situation 
the  gradient  field  of  the  coils  represents  another  independent  noise  source 
which  degrades  effectiveness  of  the  gradient  cancellation  technique. 

In  addition  to  this  the  gradient  field  of  the  coils  fluctuates  in  response 
to  variations  in  electric  currents  flowing  in  the  coils.  If  the  coil  dipole 
lies  exactly  on  the  xj  axis  of  the  gradiometer,  the  gradiometer  is  completely 
insensitive  to  dipole  orientation  so  that  current  fluctuations  in  the  coils 
produce  zero  gradiometer  response.  However,  for  a  non-perfect  alignment  coil 
current  fluctuations  produce  a  gradiometer  response  given  approximately  by  the 
product  of  the  gradient  fluctuation  and  the  angular  misalignment.  Hence  long 
term  stability  limitations  on  the  coil  power  supplies  and  even  variations  in 
the  ambient  temperature  around  the  coils  represent  other  independent  sources 
of  gradiometer  noise  which  further  degrade  effectiveness  of  gradient  cancel- 
1 ation. 

In  fact  limitations  on  the  coil  current  stability  eventually  eliminated 
the  use  of  the  gradient  cancellation  technique.  Our  standard  off-the-shelf 
power  supplies  for  the  coils  provide  a  long  term  stability  of  order  .15*  which. 
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Figure  10.  Gradiometer  noise  spectra  comparing  noise  levels  with  coils  on 

( - )  and  coils  off  ( — ).  The  coils  were  aligned  to 

cancel  tower  gradients  and  gradiometer  was  aligned  so  coil  dipole 
ifi,  was  along  x(  axis  of^g radiometer  and  gradiometer  pick-up  loops 
were  in  plane  of  ifi  and  x{.  Solid  line  is  gradiometer  noise  at  ex 
perimentally  measured  minimum  noise  position. 


for  a  steady  gradient  of  100  nT/m,  produces  gradient  fluctuations  of  150  pT/m. 
This  does  not  represent  a  significant  noise  source  when  the  coil  dipole  is 
positioned  to  lie  on  the  x^  gradiometer  axis  to  within  about  0.1  degress  and 
we  easily  achieve  this  alignment  by  positioning  the  coils  to  obtain  a  null 
gradiometer  response  when  a  large  gradient  field  from  the  coils  is  turned  on 
and  off.  Using  our  original  criteria  for  the  MNO  of  the  tower  and  aligning 
the  coil  dipole  to  cancel  the  steady  tower  gradients,  our  technique  would 
simultaneously  eliminate  gradient  noise  from  coil  current  fluctuations  and 
align  the  gradiometer  to  the  correct  MNO.  Our  later  experimental  MNO  measure¬ 
ments  however  showed  that  we  could  achieve  the  requisite  noise  suppression  simply 
by  proper  orientation  of  the  gradiometer  but  they  also  showed  that  the  true 

/v 

MNO  did  not  place  the  equivalent  dipole  on  the  x-j  gradiometer  axis.  Hence  the 
gradiometer  could  not  be  aligned  to  the  MNO  with  respect  to  the  tower  gradi- 

A 

ents  while  simultaneously  placing  the  coil  dipole  on  the  x'  gradiometer  axis. 
Consequently,  after  obtaining  the  results  shown  in  Figures  9  and  10,  we  aban¬ 
doned  subsequent  attempts  to  suppress  the  gradiometer  noise  using  the  nested 
coils  and  further  work  with  the  coils  was  limited  to  cancellation  of  the  steady 
tower  gradients  for  purposes  of  balancing  the  gradiometer  in  the  presence  of 
a  gradient  field. 


VII.  GRAD IOMETER  BALANCING 


After  acquiring  the  instrument  in  January  1976  we  began  to  perform 
field  tests  on  the  instrument  at  the  La  Posta  Astrogeophysical  Observatory. 
The  results  of  the  preliminary  tests  through  approximately  January  1977 
indicating  instrument  noise  levels  in  a  magnetically  quiet  environment  and 
describing  our  preliminary  efforts  to  balance  the  instrument  were  reported 
previously  (Podney  and  Gillespie,  1977).  After  acquiring  the  fiberglass 
gimbal ,  experiments  during  July  and  August  of  1977  were  primarily  con¬ 
cerned  with  efforts  to  balance  the  instrument  in  a  large  gradient  field 

preparatory  to  the  scheduled  tower  experiment  and  a  detailed  investigation 
of  the  balance  effects  was  delayed.  Prior  to  transporting  the  gradio- 
meter  to  the  NOSC  tower  in  September  1977,  the  instrument  was  adjusted  to 
a  balance  of  about  1  part  in  10®  at  the  La  Posta  site,  from  which  time 
the  balance  was  not  adjusted  until  after  the  completion  of  the  tower  ex¬ 
periment.  After  returning  to  La  Posta  the  balance  was  checked  and  found 
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to  be  about  2  parts  in  10  which  represents  degradation  in  the  balance 
of  more  than  a  factor  of  10  from  its  initial  state. 

Following  this  discovery  we  undertook  a  more  detailed  investigation 
of  the  balancing  characteristics  of  the  instrument.  In  this  section  we 
will  describe  quantitatively  and  in  some  detail  more  recent  efforts  to 
balance  the  instrument  using  the  new  gimbal  which  allows  more  general 
rotations  of  the  gradiometer  than  were  used  in  the  earlier  work.  We  also 
describe  some  anomolous  features  of  the  instrument  behavior  that  appear 
to  severely  limit  the  accuracy  with  which  the  instrument  can  be  balanced. 


A .  Balancing  Procedure 

All  of  our  balancing  efforts  depended  on  the  technique  of  producing  a 
null  response  from  the  gradiometer  during  a  rotation  in  a  uniform  magnetic 
field.  The  capability  to  rotate  the  instrument  about  an  arbitrary  axis  requires 
a  gimbal  mounting  having  3  independent  axes  of  rotation  as  shown  previously 
in  Figure  3.  A  0  rotation  allows  the  instrument  to  be  tilted  away  from  the 
vertical  position  through  an  angle  0,  a  $  rotation  allows  the  instrument  to 
be  rotated  360°  about  a  vertical  axis,  and  an  ft  rotation  provides  for  the 
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instrument  to  be  rotated  360°  about  its  own  xj  axis.  The  wooden  gimbal*  with 
which  we  performed  our  preliminary  experiments  was  restricted  to  0  and 
rotations  only,  but  the  fiberglass  gimbal**  which  we  received  in  July,  1977 
provided  us  with  the  additional  capability  to  perform  ft  rotations.  Both  gim¬ 
bals  were  restricted  to  a  maximum  tilt  angle  of  0  =  ±45°  to  avoid  problems  with 
the  cryogenic  dewar  and  helium  bath. 

The  balance  adjustment  on  our  instrument  consists  of  three  thin  mutually 
perpendicular  superconducting  niobium  discs  mounted  on  separate  slides  near 
the  top  pick-up  loop.  The  balance  is  then  adjusted  by  changing  position  of 
the  discs  using  rods  leading  to  the  top  of  the  dewar,  which,  in  effect,  alters 
size  and  orientation  of  the  pick-up  loops.  If  we  represent  the  imbalance  of  the 
gradiometer  by  a  vector  quantity  as  in  equation  (4),  then  in  the  gradiometer 

The  wooden  gimbal  was  built  by  Arlington  Woodworking  Company,  Arlington,  VA. 

The  fiberglass  gimbal  was  designed  by  Mechanics  Research  Inc.  of  Los  Ange¬ 
les,  CA  and  built  by  the  Naval  Ocean  Systems  Center,  San  Diego,  CA. 
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basis  each  of  the  niobium  discs  adjusts  one  component  of  the  imbalance 
vector,  t. 

Now,  in  a  uniform  field,  Sq,  which  makes  an  angle,  a,  with  a  vertical 
axis  (the  earth's  field,  for  example),  an  U  rotation  with  0=0°  produces 
a  gradiometer  response  due  to  only  the  ^  and  ^  components  of  the  imbal¬ 
ance,  since  the  component  does  not  change  its  orientation  (and  hence 

is  constant).  The  simplest  and  most  widely  used  method  for  balancing 
a  gradiometer  now  proceeds  as  follows.  The  and  ^  components  are  adjusted 
until  an  fl  rotation  with  6=0°  produces  a  null  response  from  the  grad¬ 
iometer,  a  condition  which  should  insure  that  ^  ~  0  since  the  contri¬ 

butions  from  ^  and  will  be  90°  out  of  phase  so  that  they  cannot  cancel. 
The  gradiometer  is  then  tilted  away  from  vertical  to  some  angle,  6,  and 
the  component  is  adjusted  to  give  a  null  response  during  a  360°  ro¬ 
tation  about  a  vertical  axis  (a  <p  rotation).  The  procedure  can  then  be 
iteratively  repeated  to  eliminate  effects  of  mutual  interaction  between  bal¬ 
ance  discs. 

This  is  the  technique  employed  universally  by  both  users  and  sup¬ 
pliers  of  superconducting  gradiometers.  During  our  preliminary  experi¬ 
ments  we  too  were  limited  to  this  technique  of  balancing  due  to  the  re¬ 
strictions  of  the  wooden  gimbal,  but  our  fiberglass  gimbal  provided  the  cap¬ 
ability  to  rotate  the  instrument  about  an  arbitrary  axis  and  confirm  (or  re¬ 
fute)  the  validity  of  this  definition  of  the  balance. 

B.  Balancing  Experiments  Using  3  Axes  of  Rotation 

In  January  of  1978  we  began  to  examine  balancing  characteristics 
of  the  device  more  closely.  As  an  initial  approach  to  the  problem  we  began 


by  using  the  procedure  described  above  of  adjusting  the  and  ^  components 
of  the  imbalance  vector  to  produce  a  null  response  during  an  ft  notation  at 
0=0.  We  then  tipped  the  gradiometer  away  from  vertical  to  0  =  10°  and  per¬ 
formed  another  ft  rotation,  and  finally,  to  check  the  results,  we  returned  the 
gradiometer  to  its  original  vertical  orientation  and  performed  a  third  ft  ro¬ 
tation.  The  results  are  shown  in  Figure  11  in  the  form  of  chart  recordings 
where  the  gradiometer  output  in  nT/m  is  plotted  on  the  vertical  scale  as  the 
gradiometer  is  rotated  stepwise  in  60°  increments.  The  numbered  flat  parts 
of  the  traces  represent  the  dwell  points  and  the  corresponding  values  of 
ft;  the  inset  shows  the  gradiometer  orientation  with  respect  to  the  earth  and 
the  earth's  magnetic  field  at  0  =  0°  and  0  =  10°. 

Using  the  known  magnitude  and  direction  of  the  earth's  field  we  can  de¬ 
duce  values  for  6-j  and  6,,  from  each  trace  by  finding  sin  ft  and  cos  ft  compon¬ 
ents  respectively.  (The  component  makes  only  a  dc  contribution  during  any 

given  ft  rotation.)  Fourier  decompositions  for  all  three  traces  are  given  in 
Table  I  along  with  corresponding  values  of  6-j  and  From  Table  I  it  appears 
that  both  £-|  and  components  of  balance  change  dramatically  as  a  function 
of  0  but  in  a  reproducible,  non  hysteretic  fashion  as  shown  by  the  very  simi¬ 
lar  results  for  traces  A  and  C,  both  at  0  =  0.  The  dashed  line  in  Trace  B 
shows  the  expected  gradiometer  output  for  the  values  of  6-j  and  6^  measured 
from  Trace  A.  This  dashed  line  differs  from  Trace  A  due  to  slightly  different 
gradiometer  orientation  with  respect  to  the  earth's  field  at  0  =  10°. 

To  verify  our  observation  we  modified  the  procedure  and  first  adjusted 
the  £.|  and  balance  components  to  give  a  null  response  during  an  ft  rotation 
with  0  =  15°.  Keeping  0  =  15°,  we  then  adjusted  the  component  to  minimize 
the  gradiometer  response  during  a  <J>  rotation.  In  principle  this  should  complete 
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Figure  11.  Chart  recorder  traces  showing  apparent  change  in  gradiometer  bal¬ 
ance  as  0  is  changed.  Traces  were  recorded  in  sequence  of  A,  B, 
C.  Numbered  flat  regions  mark  values  of  il  at  dwell  points  where 
gradiometer  output  was  recorded.  Spikes  arise  from  eddy  currents 
in  mylar  insulation  during  instrument  rotation.  Table  I  shows 
numerical  values  of  apparent  balance. 
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Table  1.  Fourier  decomposition  o(  ^radiometer  response  during  W  rotations 
and  corresponding  values  of  <$i  and  deduced  from  sin  V  and 

cos  W  Fourier  components. 
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the  balancing  process  once  enough  iterations  are  performed  to  eliminate  inter¬ 
actions  between  balance  disks.  We  then  performed  two  ft  rotations;  the  first 
with  9  =  15°  the  second  with  9=0°  with  the  result  shown  in  Figure  12.  Again 
the  Fourier  decompositions  and  corresponding  values  of  <5-|  and  62  are  shown 
in  Table  I  and  the  dashed  line  in  Trace  B  shows  the  expected  response  on  the 
basis  of  the  6^  and  §2  determined  from  Trace  A.  Thus  after  carefully  balanc¬ 
ing  the  gradiometer  to  1  part  in  106  with  respect  to  all  three  axes,  we  see 
that  the  apparent  and  ^  components  of  the  balance  change  dramatically 
when  the  dewar  is  tilted  to  a  different  value  of  9.  The  change  in  balance 
shown  in  Figure  12  represents  a  change  in  +  ^1  ^rom  a  value  at  9  =  15° 

of  about  2x10  ®  m"^  to  a  value  of  about  1x10"^  m’"*  at  9  =  0°.  In  short  it 
seems  clear  that  while  it  is  relatively  straightforward  to  balance  our  grad¬ 
iometer  to  a  part  in  10^  using  the  common  definition  of  the  term  "balance", 
the  balance  achieved  using  a  two-axis  rotation  scheme  does  not  necessarily 
hold  for  an  arbitrary  rotation  axis. 

The  source  of  this  behavior  is  still  unknown  although  we  suspect  that 
the  effect  may  be  the  result  of  effective  changes  in  the  level  of  the  helium 
bath  during  the  rotation  about  a  non- vertical  axis.  Since  the  helium  itself 
is  diamagnetic  and  the  G-10  fiberglass  in  the  dewar  is  thought  to  contain  para¬ 
magnetic  impurities,  rotations  of  the  gradiometer  about  a  non-vertical  axis 
or  changes  in  the  tilt  angle,  0,  can  bring  different  portions  of  the  dewar  into 
contact  with  the  helium  bath  producing  temperature  and  magnetic  fluctuations 
in  the  dewar  walls  as  well  as  magnetic  fluctuations  in  the  gradiometer  due  to 
changes  in  gradiometer  orientation  with  respect  to  the  surface  of  the  helium 
bath.  These  effects  could  lead  to  apparent  non-hysteretic  changes  in  the 
balance  for  ft  rotations  at  different  values  of  9  as  shown  in  Figure  11.  Since 
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Figure  12.  Chart  recorder  traces  showing  apparent  change  in  balance  as  0  is 
changed.  All  3  axes  of  gradiometer  were  initially  balanced  at 
6  =  15°,  then  we  recorded  an  ft  rotation  at  0  -  15°  (trace  A)  and 
an  ft  rotation  at  0  =  0  (trace  B).  Change  in  balance  is  shown 
numerically  in  Table  I.  Numbered  flat  regions  mark  values  of  ft 
at  dwell  points. 
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our  measurements  used  a  very  rigid  stationary  platform,  this  did  not  repre¬ 
sent  a  limitation  on  our  experiment.  For  applications  requiring  measurements 
from  moving  platforms  however,  such  anamolous  balance  behavior  could  seriously 
impair  the  capabilities  of  the  instrument.  Although  the  effect  could  be  pe¬ 
culiar  to  our  particular  device,  further  work  with  other  instruments  using  3 
independent  axes  of  rotation  is  highly  desirable. 

We  should  also  point  out  that  the  behavior  depicted  in  Figures  11  and 
12  cannot  be  explained  by  the  presence  of  a  local  anamolous  gradient  in  the 
earth's  field.  The  argument,  which  involves  an  additional  experimental  ob¬ 
servation,  is  presented  in  Appendix  A. 

C.  Hysteretic  Effects  in  Balancing 

In  addition  to  the  effect  described  above  we  also  found  a  very  large  hy¬ 
steretic  effect  which  first  appeared  during  January,  1978.  The  problem  is  char¬ 
acterized  by  large  and  very  sudden  shifts  in  the  dc  output  voltage  of  the  grad- 
iometer  as  the  device  is  tilted  to  large  tilt  angles.  This  behavior  is  ac¬ 
companied  by  large  changes  in  the  apparent  ^  and  balance  components  when 
the  gradiometer  is  returned  to  a  vertical  orientation.  The  behavior  of  the 
instrument  during  0  rotations  is  shown  in  Figure  13  again  in  the  form  of  chart 
recordings  where  the  numbered  flat  parts  of  the  trace  represent  the  values  of 
0  at  the  dwell  points  during  a  stepwise  rotation.  The  arrows  labeled  A,  C, 
and  D  mark  the  large  dc  shifts  observed  as  the  dewar  was  tilted  to  angles  larger 
than  30°.  The  arrow  labeled  B  marks  a  smaller  shift  on  going  from  15°  to  30° 
in  0.  (Note  that  this  effect  is  not  the  same  as  the  hysteretic  effect  described 
by  Podney  and  Gillespie,  1977  which  was  observed  even  during  <J>  rotations  with 
0=0°  and  was  much  smaller.)  Figure  14  shows  the  apparent  effect  of  these 
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Figure  14.  Chart  recordings  showing  apparent  change  in  balance  as  observed  at 
0=0°  after  tilting  dewar  30°  from  vertical.  Traces  were  recorded 
in  sequence  A,  B,  C.  Numbered  flat  regions  mark  values  of  ft  at 
dwell  points  where  gradiometer  output  was  recorded. 
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dc  shifts  on  the  and  <$,,  components  of  the  balance  vector  of  the  instrument. 
Trace  A  shows  the  instrument  response  for  an  ft  rotation  with  0  =  0°,  trace  B 
is  the  response  for  an  ft  rotation  but  with  0  =  30°  and  trace  C  is  from  a  second 
ft  rotation  at  9  =  0°.  The  Fourier  decomposition  of  the  three  traces  is  given 

A  A 

in  Table  I  showing  that  the  component  of  6  in  the  (xj.x^)  plane  has  increased 
by  a  factor  of  about  2. 

Considerable  time  and  effort  were  expended  in  attempting  to  understand 
this  effect,  and  although  its  origin  is  still  unknown,  the  large  sudden  shifts 
suggest  some  type  of  a  motion  which  occurs  at  large  tilt  angles.  During  the 
period  in  which  the  gradiometer  was  warmed  up  we  did  find  that  a  .8  mm  diameter 
fiberglass  rod  which  maintains  tension  on  the  balance  slides  to  prevent  slip- 

A 

page  had  cracked  partly  releasing  tension  on  the  x^  balance  slide.  However, 
after  the  fiberglass  rod  was  replaced  and  the  instrument  again  cooled  down 
to  4°K,  we  still  observed  the  sudden  hysteretic  dc  shifts.  Furthermore,  we 
confirmed  during  the  subsequent  warmups  required  to  readjust  the  permanent 
balance  disks  that  the  new  fiberglass  rod  was  still  intact  and  providing  ade¬ 
quate  tension  on  all  three  of  the  balance  slides. 

An  alternative  explanation  has  been  suggested*  concerning  the  thermal 
equalization  shield  in  the  vacuum  jacket  of  the  dewar.  Thermal  equali¬ 
zation  along  the  vertical  axis  of  the  dewar  is  achieved  with  several  layers 
of  copper  "coil  foil"  constructed  of  layers  of  copper  wire  held  together 
with  some  type  of  glue  or  adhesive  varnish.  The  "coil  foil"  is  then  im- 


*  This  explanation  was  proposed  during  a  private  conversation  with  another 
commercial  superconducting  instrument  supplier  who  is  familiar  with  the 
dewar  design  used  by  SCT. 


bedded  between  layers  of  superinsulation.  If  no  rigid  mechanical  support 
is  provided  for  this  assembly  in  the  vacuum  jacket  of  the  dewar  (which 
we  believe  to  be  the  case)  then  tilting  the  dewar  through  large  angles 
could  easily  cause  the  "coil  foil"  and  superinsulation  wrapping  to  crush 
down  against  one  side  or  the  other  of  the  inner  wall  of  the  vacuum  space. 
Unfortunately,  this  hypothesis  can  be  verified  only  by  the  costly  process 
of  completely  disassembling  the  dewar,  rebuilding  the  interior  of  the  dewar's 
vacuum  jacket,  and  conducting  further  tests. 

It  should  be  emphasized  here  that  the  hysteretic  effect  just  described 
was  observed  only  for  tilt  angles  greater  than  about  25°  and  is  not  the 
same  as  the  effect  displayed  in  Figure  11  which  was  distinctly  not  hysteretic. 
We  believe  that  the  large  hysteretic  effect  is  probably  due  to  some  mechan¬ 
ical  motion  peculiar  to  this  particular  instrument  and  could  be  avoided  with 
improved  design.  The  definition  of  the  balance  with  respect  to  two  or  three 
axes  of  rotation  however,  represents  a  more  fundamental  problem  since  the  two- 
axis  system  is  widely  accepted  as  defining  the  correct  level  of  balance.  Mag¬ 
netic  effects  from  the  dewar  walls  and  the  helium  bath  can  be  minimized  by 
proper  instrument  design  but  the  lack  of  a  definitive  answer  to  the  question 
of  instrument  balance  for  arbitrary  rotations  represents  a  serious  obstacle 
to  the  use  of  this  type  of  instrument  on  nonstationary  platforms. 

As  a  consequence  of  the  effects  described  above, it  was  difficult  to  de¬ 
fine  a  unique  level  of  balance  for  the  instrument.  In  practice,  for  the  1978 
oceanographic  experiments  we  balanced  the  instrument  to  a  level  of  about  1  part 
in  10^  using  the  conventional  technique  which  then  produced  an  apparent  balance 

C 

for  an  ft  rotation  of  about  5  parts  in  10  for  a  tilt  angle  of  15°.  We  then 
moved  the  instrument  to  the  NOSC  tower  after  providing  a  mechanical  means  of 
restoring  the  balance  to  a  level  of  about  1  part  in  105  in  the  event  of  a  large 


change  in  balance  occurring  during  transportation.  The  mechanical  arrange¬ 
ment  consisted  of  small  brass  collars  attached  to  the  balancing  rods  at  the 
top  of  the  dewar  which  allowed  us  to  return  each  of  the  niobium  balance  discs 
to  a  certain  mechanical  position  with  respect  to  the  gradiometer  pick-up  loop. 
The  desired  position  was  initially  determined  by  balancing  the  gradiometer 
after  which  the  brass  collars  were  attached  such  that  the  position  of  the  nio¬ 
bium  discs  could  be  restored  if  needed.  For  our  modest  requirements  on  the 
balance  level  this  approach  proved  adequate. 
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VIII.  GRAD IOMETER  ALIGNMENT  AND  BALANCING  IN  A  GRADIENT  FIELD 


The  motivation  for  some  of  the  work  discussed  in  this  section  was 
provided  by  our  desire  to  have  the  capability  to  balance  the  gradiometer 
in  the  presence  of  a  large  gradient,  as  well  as  to  suppress  the  noise  from 
the  ambient  gradient  field.  The  mechanical  nature  of  the  balancing  mech¬ 
anism  raised  questions  concerning  the  stability  of  the  balance  during 
transportation  of  the  gradiometer  from  the  La  Posta  site  to  the  NOSC  tower, 
but  without  some  means  to  separate  the  gradient  and  imbalance  components 
of  the  gradiometer  response  in  the  presence  of  the  gradient  field  of  the 
tower  there  was  no  way  to  adjust  or  even  check  the  instrument  balance 
once  it  was  installed  on  the  tower.  However,  if  in  some  manner,  the  grad¬ 
ient  field  of  the  tower  could  be  precisely  determined,  then  a  coordinate 
system  coincident  with  the  principal  axes  of  that  gradient  field  would 
provide  a  reference  frame  in  which  the  gradient  response  would  be  zero 
for  selected  rotations  of  the  gradiometer.  Alternatively,  the  tower  grad¬ 
ients  might  be  cancelled  using  the  nested  coils  to  allow  the  instrument 
to  be  balanced  by  rotations  about  arbitrarily  selected  axes.  And  finally, 
of  course,  a  precise  knowledge  of  the  steady  gradient  field  would  allow 
a  similarly  precise  alignment  of  the  gradiometer  to  suppress  noise  gener¬ 
ated  by  relative  motions  of  the  gradiometer  and  the  ambient  steady  grad¬ 
ients. 

It  is  in  principle,  possible  to  completely  separate  the  gradient 
and  imbalance  contributions  to  the  gradiometer  response  during  a  series 
of  rotations  of  the  gradiometer  in  the  presence  of  a  magnetic  field  having 


spatial  gradients  as  well  as  a  spatially  uniform  component.  As  discussed 
later  in  this  section,  the  value  of  the  balance  determined  for  our  instru¬ 
ment  by  this  calculation  was  unreliable  and  in  any  event  the  calculation 
was  not  considered  sufficiently  precise  to  allow  balance  adjustments  to 
the  requisite  level.  Nonetheless  a  precise  measurement  of  the  gradient 
field  of  the  tower  should  provide,  to  the  same  precision,  the  position, 
orientation,  and  strength  of  the  equivalent  dipole  which  would  allow  us 
to  then  cancel  the  steady  gradients  using  the  nested  coils.  A  residual 
gradient  strength  after  cancellation  of  order  2  nT/m,  should  subsequently 
permit  balancing  to  a  few  parts  in  10*\  To  achieve  this  level  of  cancel¬ 
lation  the  polar  and  azimuthal  angles  giving  the  position  and  orientation 
of  the  equivalent  dipole  must  be  measured  to  an  accuracy  of  about  0.5°. 

The  concept  of  precise  alignment  of  the  gradiometer  with  a 
gradient  field  to  suppress  noise  in  the  gradiometer  due  to  the  presence 
of  the  steady  gradients  was  introduced  in  Section  V.  To  eliminate  the  noise 
generated  by  effective  rotations  of  the  gradiometer  due  to  twisting  or 
oscillations  of  the  boom  we  built  a  set  of  nested  coils  with  which  we  hoped 
to  cancel  the  gradient  field  of  the  tower  and  provide  a  gradient  free  re¬ 
gion  of  space  at  the  site  of  the  gradiometer.  To  suppress  noise  generated 
by  small  perturbations  in  the  gradient  field  of  the  tower  we  proposed 
aligning  the  gradiometer  such  that,  to  first  order,  the  gradiometer  response 
to  small  fluctuations  of  the  gradient  field  would  be  minimized.  Both  of 
these  techniques  were  proposed  in  the  context  of  representing  the  gradient 
field  of  the  tower  in  the  form  of  its  equivalent  dipole  (as  described  in 
section  IV)  and  using  the  position  and  orientation  of  the  equivalent 
dipole  to  align  the  gradiometer,  the  coils,  and  the  coil  dipole  moment. 
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Since  the  noise  introduced  into  the  gradiometer  is  to  first  order  propor¬ 
tional  to  its  angular  misalignment  in  the  gradient  field,  the  effective¬ 
ness  of  both  techniques  depends  upon  the  ability  to  determine  precisely 
the  parameters  of  the  equivalent  dipole  (or  alternatively  the  five  inde¬ 
pendent  elements  of  the  gradient  matrix). 

We  initially  envisioned  the  alignment  procedure  as  consisting  of 
two  steps;  the  first  to  measure  the  approximate  parameters  of  the  equiva¬ 
lent  dipole  (position,  orientation,  and  strength)  and  the  second  to  define 
their  values  with  greater  precision  to  allow  the  final  alignment  adjust¬ 
ments.  The  first  estimate  of  the  equivalent  dipole  parameters  was  made 
by  essentially  "mapping  out"  the  ambient  magnetic  and  gradient  fields 
using  a  systematic  set  of  gradiometer  rotations  then  fitting  the  gradio¬ 
meter  output  to  the  corresponding  mathematical  model.  We  expected  this  method 
to  provide  alignment  values  to  within  about  5°  of  the  true  value.  To 
provide  the  additional  required  precision  to  bring  the  gradiometer  and 
coil  alignment  to  within  about  1°  we  initially  proposed  a  method  by  which 
we  hoped  to  locate  the  principle  axes  by  observing  the  gradiometer  output 
during  a  series  of  0  rotations.  In  practice  this  process  failed  and  we 
developed  alternative  approaches  to  the  alignment  problem,  as  we  will 
describe  in  the  following  paragraphs. 

A .  Initial  Estimate  of  A  Gradient  Field 

To  provide  a  first  estimate  of  the  ambient  gradients  we  used  a  tech¬ 
nique  based  on  our  model  for  the  gradiometer  response  under  rotations 
in  a  gradient  field  as  given  by  equation  (15).  The  actual  measurements 
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typically  consisted  of  360°  stepwise  rotations  in  4>  and  ft  using  60°  and 
90°  steps  respectively  and  a  0  rotation  from  -20°  to  +20°  in  steps  of 
5°.*  In  effect,  at  each  value  of  0  and  <{>  we  perform  a  stepwise  ft  rota¬ 
tion  through  360°  recording  the  gradiometer  output  at  each  step,  a  process 
which  provides  a  pointwise  mapping  over  the  entire  range  of  0,  $,  and  ft. 
The  calculation  of  the  gradient  field  and  imbalance  vector  is  then  per¬ 
formed  as  follows.  First  we  write  the  gradiometer  response,  from  equa¬ 
tion  (15)  as 

r(0, 4>,ft)  =  C(0,4>)  +  A(0,<j>)cos  ft  +  B(0,<{i)sin  ft  .  (34) 

We  can  now  write 

A(0,4>)  =  A  (0)  +  A^Ojcos  <})  +  A2(9)sin  <J> 

+  A3(0)cos  2<J>  +  A4(0)sin  2<J>  (35) 

and  similarly  for  B(G)  and  C(0),  and  finally  each  of  the  coefficients 
C^(0-|)  A^ ( 0 )  and  B..(0)  can  Fourier  decomposed  as,  for  example 

A^(0)  =  A?  +  a!  cos  0  +  A?  sin  0  +  A?  cos  20  +  A^  sin  20  (36) 

The  explicit  form  for  each  of  the  constants  can  now  be  found  simply  by 
organizing  the  expression  in  equation  (15)  as  Fourier  expansions  in  the 
above  manner.  The  coefficients  A'j,  B?,  C'?  are  now  functions  only  of  the 
five  independent  elements  of  the  gradient  matrix,  the  g.  .'s  and  the  three 

■  J 

Although  the  gimbal  allows  a  range  of  -45°  <_  0  <_  45°,  the  dynamic  range 

of  the  gradiometer  is  insufficient  to  allow  rotations  at  tilt  angles 

greater  than  20°  in  the  large  gradient  field  of  the  tower.  Limitations 

due  to  the  balance  instability  dictated  a  similar  restriction  for  our 

instrument  even  in  the  absence  of  gradient  fields. 
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components  of  both  the  uniform  earth’s  field,  and  the  imbalance  vec¬ 
tor,  i.  Here  again  the  g^'s  and  the  three  components  of  bQ  are  specified 
with  respect  to  the  earth  basis  while  the  components  of  t  are  regarded 
as  fixed  with  respect  to  the  gradiometer.  Once  the  constants  are  known 
and  given  the  calibration  constant  for  the  gradiometer  and  the  magnitude 
of  the  earth's  field,  it  is  trivial  to  calculate  the  components  of  the 
gradient  field,  the  imbalance,  and  even  the  individual  components  of  t5Q. 

The  calculational  approach  is  now  straightforward.  First  the  data  are 
organized  into  individual  fl  rotations,  each  of  which  is  decomposed  into 
its  Fourier  components  to  give  values  for  A(0,<f>),  B(0,<J>),  and  C(0,<f>). 

In  exactly  the  same  manner  a  second  Fourier  decomposition  in  <j>  gives  the 

O 

coefficients  C^(o),  A^(o),  and  IL (o)  and,  in  turn,  a  Fourier  analysis  in 
6  gives  the  A^'s,  B'j's,  and  C^'s.  The  computer  program  which  performs 
these  calculations,  called  GRADB,  is  described  in  Appendix  B,  as  is  a  brief 

O 

description  of  the  tests  we  performed  to  determine  the  consistency  of  the 
calculation  and  its  sensitivity  to  the  presence  of  random  noise  in  the 
data  collected  from  the  gradiometer. 

The  GRADB  program  was  used  extensively  both  at  La  Posta  and  during 
all  of  the  tower  experiments  to  provide  a  first  estimate  of  the  ambient 
gradient  field.  At  the  La  Posta  site  we  performed  gradiometer  rotations 
in  an  essentially  gradient  free  environment  (the  earth's  field)  as  well 
as  in  the  gradient  fields  of  the  nested  coils  and  the  iron  sphere  used 
to  simulate  the  tower  gradients.  The  results  of  the  experiments  at  La 
Posta  are  shown  in  Table  II,  and  those  from  the  tower  in  Table  III.  The 
rotations  performed  in  the  earth's  field  alone  gave  quite  reproducible 
••  .lilts  for  the  gradients  and  the  imbalance  vector  but  the  numbers  calcu- 
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Table  II. 


Numerical  results  of  GRADB  calculation  using  data  collected 
during  various  sets  of  systematic  gradiometer  rotations  at  the 
La  Posta  site. 


Date 

Gradient  Field  of 

g 

gll  /g12/913/g22/g32 

VV63 

Bi/b2/b. 

(NT/m) 

no"  V) 

(NT) 

13  Jan  78 

earth 

-.592 

3.37 

17441 

-.291 

-9.71 

19218 

-.028 

-1.355 

-7.48 

36762 

-.1436 

17  Jan  78 

earth 

-1.500 

4.06 

16896 

.278 

-15.1 

12302 

-.339 

-2.185 

-3.69 

39852 

-.009 

19  Jan  78 

earth 

2.2 

-1.13 

4.41 

18314 

-.280 

-11.74 

13571 

+  .002 
-1.90 

-4.53 

38800 

.161 

18  Jan  78 

coils 

26.3 

-.170 

13.12 

-29798 

Mx=° 

-.297 

-18.33 

26390 

M  =0 

24.34 

-6.95 

20991 

y 

M  =25  nT/m 

-3.02 

z 

-22.33 

18  Jan  78 

coi  1  s 

22.7 

20.19 

9.09 

22267 

Mx=° 

4.12 

-12.25 

16816 

M  =25  nT/m 

.593 

-2.78 

35304 

M  =0 

-44.63 

12  April  78 

z 

iron  sphere  at 
IT 

10.08 

171.6 

-29.34 

112.98 

-31641 

-40.04 

-118.49 

30323 

+25.31 

+73.11 

6.79 

10217 

Table  II.  (cont.) 


Date  Gradient  Field  of  9  9n  /9]2/9VJg22/93Z  V62/63  Bl/B2/B3 


23  April  78 

iron  sphere 

22.7 

20.54 

3.83 

-3728 

at  12  ft  with 

-.249 

56.12 

10342 

coils  at 

3.33 

6.38 

43636 

Mx=27.9  nT/m 

14.83 

My=50.0  nT/m 

-5.52 

Mz=54.5  nT/m 

25  April  78 

iron  sphere 

26.5 

13.45 

15.14 

6089 

at  12  ft. 

-.315 

108.72 

23076 

coils: 

3.45 

-3.91 

-391077 

Mx=29.2  nT/m 

5.08 

My =4 3. 2  nT/m 

-27.54 

Mz=33.4  nT/m 

26  April  78 

iron  sphere 

116.7 

-16.28 

43.29 

-7952 

at  12' 

-27.66 

192.4 

11215 

8.49 

209.50 

42849 

89.44 

-83.29  D 


.1 


f 
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Table  III. 


Numerical  results  of  GRAOB  calculation  using  data  collected 
during  gradiometer  rotations  in  the  steady  gradient  field  of 
the  NOSC  tower 


Date  gradient  field  of 


16  May  78  NOSC  Tower 


18  May  78  NOSC  Tower 


19  May  78  NOSC  Tower 


31  May  78  NOSC  Tower 


5  June  78  NOSC  Tower  with 

coils  at 
Mx=19.9  NT/m 

My=48. 9  NT/m 

Mz=l 25.7  NT/m 


g 

(nT/m) 


911  /912/g13/g22/g23 


238.8  -123.6 

-28.6 
15.8 

332.5 
-41.1 

224.0  -135.2 

-28.5 
10.2 
333.7 
-41.6 

240.7  -130.9 

-28.8 

8.8 

343.5 
-42.3 


208.1  -141.0 

-0.6 
45.2 
327.5 
-43.2 


6.9 

-3.7 

0.1 


-3.6 


lated  for  the  steady  field  of  the  earth  are  definitely  incorrect.  For  our 
initial  gradiometer  orientation  with  the  earth  and  gradiometer  bases  coinci¬ 
dent  we  would  expect  to  be  near  zero  and  =  2B£.  Although  the  final  re¬ 
sult  of  the  calculation  is  wrong,  its  reproducibility  suggests  the  presence 
of  a  systematic  error  in  either  the  data  or  the  calculation.  Since  the  tests 
we  performed  on  the  GRADE  program  effectively  eliminated  the  calculation  as 
a  possible  source  we  concluded  that  the  gradiometer  response  under  our  sys¬ 
tem  of  rotations  contained  a  systematic  deviation  from  the  mathematical  model 
we  were  using.  Subsequent  experiments  on  the  balancing  characteristics  of  the 
gradiometer,  which  we  described  in  the  last  section,  confirmed  that  the  behav¬ 
ior  of  the  gradiometer  imbalance  was  much  more  complicated  than  originally 
anticipated. 

Measurements  in  a  known  gradient  field  are  shown  in  the  data  from  18 
January  for  a  set  of  rotations  in  the  gradient  field  of  the  nested  coils. 

For  both  cases  the  strength  of  the  equivalent  dipole,  g,  is  approximately 
correct  and  the  magnitudes  of  the  respective  gradient  elements  are  about 
as  expected;  that  is,  large  longitudinal  gradients  for  the  case  where 
My  is  nonzero  and  large  transverse  elements  for  a  nonzero  .  In  general, 
these  results  provided  a  certain  level  of  confidence  in  the  values  cal¬ 
culated  by  the  GRADB  routine  for  the  gradients  and  simultaneously  demon¬ 
strated  that  the  calculations  for  the  components  of  the  imbalance  and 
magnetic  fields  were  totally  unreliable. 

This  conclusion  may  appear  somewhat  remarkable  since  initially  the  var¬ 
ious  components  of  the  gradiometer  response  are  totally  mixed  together  into  a  sin¬ 
gle  output.  Nonetheless,  since  the  gradients  give  a  response  with  contains  2$  and  20 


components  in  the  Fourier  decomposition,  the  model  can  theoretically 
differentiate  completely  between  a  gradient  response  and  contributions 
due  to  imbalance  or  other  effects  which  produce  no  20  or  2<p  components. 

This  statement  must  be  qualified  somewhat  due  to  the  experimental  limitation 
imposed  by  the  cryogenic  nature  of  the  instrument  and  its  helium  bath 
which  restricts  the  tilt  angle  to  the  range  -45°  <_  0  _<  45°.  For  a  ro¬ 
tation  in  4> ,  the  <J>  and  2<t>  Fourier  components  are  orthogonal  over  the  range 
of  rotation,  namely  0  _<<{>_<  2it,  and  are  easily  separated  in  the  GRADB 
calculation,  but  over  the  accessable  range  of  the  0  rotation  the  0  and 
20  Fourier  components  are  not  orthogonal.  In  fact  the  functional  forms 
of  cos  0  and  cos  26  over  this  domain  are  very  similar  as  are  those  for 
sin  0  and  sin  20.  Consequently  any  systematic  0  dependence  which  is  not 
accounted  for  in  our  model  has  the  potential  for  introducing  significant 
errors  into  the  GRADB  calculation  of  both  gradients  and  imbalance.  This 
is  precisely  the  type  of  behavior  we  observed  during  our  experiments  with 
the  gradiometer  balance  described  in  Section  VII.  Nonetheless,  as  shown 
in  Table  II  our  measurements  on  18  January  in  the  known  gradient  field 
of  the  nested  coils  indicate  that  the  gradient  calculation  is  not  too 
seriously  impaired  at  least  in  a  gradient  field  of  order  25  nT/m.  In 
contrast  calculations  from  the  same  data  produced  excessively  large  values 
for  the  components  of  the  imbalance  vector  $  and  obviously  incorrect  values 
for  the  earth's  field,  (also  shown  in  Table  II)  which  are  undoubtedly 
a  reflection  of  the  anamolous  behavior  of  the  gradiometer  balance.  If  so, 
the  consistency  of  the  results  from  our  gradiometer  rotations  during 
January  indicate  that  the  effect  is  very  reproducible 


and  probably  non-hysteretic,  which  again  describes  the  observed  balance 
behavior  for  tilt  angles  less  than  about  25°.  As  a  final  comment  on  this 
discussion,  although  the  GRADB  calculation  produces  approximately  the 
correct  value  for  a  gradient  field  of  about  25  nT/m,  we  still  expect  the 
deviations  in  the  balance  to  limit  the  accuracy  of  the  gradient  calcula¬ 
tions  at  some  lower  value  of  the  gradient  field. 

We  should  also  note  here  that  our  discovery  of  the  stainless  steel 
screws  near  the  pick-up  loops  gave  rise  to  the  hope  that  subsequent  experi¬ 
ments  would  provide  better  results  for  the  calculation  of  the  gradiometer 
balance  but  experiments  conducted  after  removing  the  stainless  steel  screws 
showed  that  the  balance  still  could  not  be  described  by  any  simple  model 
consisting  of  a  simple  magnetometer- 1  ike  response  under  rotations.  Con¬ 
sequently,  we  continued  to  accept  the  gradient  calculations  as  a  rough 
estimate  of  the  true  gradient  field  and  to  disregard  the  calculations  for  the 
components  of  3  and  SQ.  The  consistency  of  the  results  obtained  on  the 
tower,  shown  in  Table  III,  provided  the  final  justification  for  this  con¬ 
clusion,  in  particular,  the  data  from  5  June  in  which  the  coils  were  aligned 
on  the  basis  of  a  GRADB  calculation  to  cancel  the  steady  tower  gradients. 

A  rough  estimate  of  the  accuracy  of  the  alignment  is  made  by  comparing 
the  gradient  strength  of  the  equivalent  dipole  from  the  tower  alone  with 
the  strength  of  the  residual  dipole  when  the  coils  are  set  to  cancel  the 
tower  gradients.  If  we  compare  the  data  from  5  June  with  an  average  value 
from  all  of  the  previous  rotations  during  May  we  find  that  the  residual 
dipole  strength  is  roughly  a  factor  of  18  smaller  than  the  uncancelled 
gradient  strength  from  the  tower.  For  a  crude  first  order  estimate  in 
which  the  residual  gradients  are  roughly  proportional  to  the  angular  mis¬ 
alignment,  this  corresponds  to  an  alignment  error  of  order  3°,  in  the 
various  angles;  about  the  expected  performance  of  the  GRADB  calculation. 


From  this  point  we  intended  to  use  a  different  technique  to  provide  a  more 
precise  alignment  with  the  gradient  field. 

B.  Precise  Alignment  with  the  Gradient  Field 

As  mentioned  briefly  in  the  introduction  to  this  section  we  initially 
proposed  a  scheme  by  which  we  hoped  to  determine  the  orientation  of  the  prin¬ 
ciple  axes  of  the  gradient  field  of  the  tower  to  an  overall  accuracy  of  about 
.5°.  This  precision  should  allow  cancellation  of  the  tower  gradients  to  a 
residual  value  of  order  2  nT/m,  which  would  subsequently  allow  alignment  and 
balancing  of  the  gradiometer  to  a  level  sufficient  to  eliminate  the  noise 
from  the  gradient  field  and  imbalance  effects.  To  perform  the  measurement 
to  the  precision  required,  we  proposed  to  start  from  the  gradients  from  the 
GRADB  calculation  and  compute  the  approximate  position  and  orientation  of  the 
equivalent  dipole  using  the  DIPOL  routine  (Appendix  C).  We  then  proposed  to 
use  360°  ft  rotations  at  small  increments  in  both  0  and  4>  and  from  equation 
(34)  determine  A (e »<f>)  and  B(0,<fi)  for  each  ft  rotation,  that  is,  at  each  value 
of  9  and  <J>.  The  0  and  <p  dependences  of  A(0,4>)  and  B(0,<|>)  can  then  be  measured 
in  the  region  near  the  principal  axes' and  compared  to  mathematical  expressions 
for  A(0,<{>)  and  B ( 0 ,<{> )  near  the  principal  axes.  (Podney  and  Gillespie,  1977, 
discuss  the  mathematical  basis  of  this  approach  in  some  detail.)  However, 
when  we  performed  the  measurements  we  found  that  the  calculation  could  not 
separate  the  imbalance  and  gradient  components  of  the  response  over  a  small 
range  of  0  and  <J>.  This  problem  is  partly  a  result  of  the  fact  that  both  imbal¬ 
ance  and  gradients  generate  only  sin  ft  and  cos  ft  responses  so  that  for  any  given 
ft  rotation  there  is  no  way  to  separate  the  imbalance  and  the  gradient  compon¬ 
ents,  which  consequently  must  be  done  on  the  basis  of  variations  in  A(0,<J>)  and 


B(e,$)  over  a  small  range  of  0  and  <j>  near  the  principal  axes.  Here  again 
the  balance  behavior  apparently  introduced  unexpected  contributions  to 
the  O  dependences  which,  when  coupled  with  the  small  range  of  0  and  <j>, 
effectively  prevented  separation  of  the  gradient  and  imbalance  components 
of  the  gradiometer  output.  As  a  result,  we  found  that  when  attempting  to 
balance  the  gradiometer  in  the  presence  of  a  large  gradient  field,  the 
procedure  could  lead  to  a  condition  in  which  a  large  gradient  component 
could  be  effectively  cancelled  by  a  similarly  large  imbalance.  Indeed, 
during  our  experiments  at  the  La  Posta  site  this  situation  occurred  and 
ultimately  proved  fatal  to  this  approach. 

After  concluding  that  this  method  would  not  yield  the  necessary  pre¬ 
cision  for  aligning  the  gradiometer,  we  developed  the  more  empirical  align¬ 
ment  technique  described  in  Section  VI  to  provide  the  requisite  noise  sup¬ 
pression,  and  further  efforts  to  precisely  measure  the  ambient  gradients 
and  corresponding  principal  axes  were  directed  toward  the  problem  of  bal¬ 
ancing  the  gradiometer  in  the  presence  of  a  large  gradient  field.  An  al¬ 
ternative  approach  to  this  problem  suggested  itself  when  we  found  that  the 
MNO  could  be  determined  to  within  about  a  degree  using  a  direct  measure¬ 
ment  of  the  gradiometer  noise  spectrum.  Once  the  MNO  had  been  measured 

A 

the  equivalent  dipole  should  lie  along  the  x-j'  axis  of  the  gradiometer 
and  the  pick-up  loops  would  define  the  plane  in  which  the  equivalent  dipole 
moment,  m,  would  lie.  This  idea  was  also  abandoned  however  when,  after 
repeated  measurements  of  the  MNO,  it  became  clear  that  the  MNO  measured 
using  the  method  described  in  Section  VI  did  not  agree  with  the  MNO  cal¬ 
culated  from  the  data  collected  during  the  gradiometer  rotations.  For 
example,  from  experiments  at  La  Posta  using  the  neste  1  coils  we  knew  that 


the  GRADB  calculations  gave  reasonable  estimates  for  the  steady  grad¬ 
ients  and  their  equivalent  dipole  yet  the  MNO  calculated  from  these  data 
disagreed  markedly  from  that  determined  by  the  direct  measurements.  This 
discrepancy,  which  will  be  discussed  in  more  detail  in  Section  IX,  was 
not  understood  until  near  the  end  of  the  August  experiment,  at  which  time 
internal  wave  measurements  performed  during  daylight  hours  prevented  ad¬ 
ditional  experiments  to  verify  our  result. 

During  the  August  internal  wave  measurements  we  made  one  further 

attempt  to  achieve  a  precise  coil  alignment  using  an  iterative  technique 
based  on  the  assumption  that  the  gradients  determined  using  the  GRADB 
calculation  are  reasonably  accurate.  If  we  perform  a  complete  set  of 
gradiometer  rotations  in  all  three  Euler  angles,  0,  <}>,  and  Q  (where  as 
before  the  tilt  angle,  0,  is  restricted  to  -20°  <  Q  ±  20°)  we  can  compute 
the  approximate  equivalent  dipole  of  the  tower  gradients  and  align  the 
nested  coils  to  cancel  it.  But  due  to  errors  in  aligning  the  coils  there 
will  be  residual  gradients  after  turning  the  coils  on.  If  the  uncancelled 
tower  gradients  are  given  by  the  gradient  matrix  Gq  and  the  residual  grad¬ 
ients  by  a  similar  matrix,  6G  then  we  can  write 

6G  =  G0  +  (-G0+6G)  (37) 

where  the  gradients,  Gc,  generated  by  the  coils  are  given  by 

Gc  =  -Go+6G  (38) 

Assuming  that  the  coil  alignment  is  accurate  to  a  few  degrees,  if  we  now 
ask  what  small  adjustments  in  the  position,  orientation,  and  strength  of 
the  coil  dipole  are  required  to  nullify  the  residual  gradients,  5G,  the 
calculation  proceeds  as  follows.  We  specify  the  position  of  an  equiva- 
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lent  dipole  of  a  gradient  field  by  the  polar  and  aximuthal  angles  <J>r  and 
6r,  its  orientation  by  the  polar  and  azimuthal  angles  4>m  and  0m>  its  strength 
by  g  and  let  this  set  of  five  parameters  be  represented  by  a  =  { <t> r , 0 r , «Pm  4>m,g}, 
Furthermore,  let  aQ  represent  the  parameters  specifying  the  equivalent 
dipole  of  the  tower  gradients,  GQ,  and  6ot  represent  a  small  arbitrary  per¬ 
turbation  of  the  parameters  a.  We  can  now  expand  the  coil  gradients,  Gc> 
in  a  Taylor  series  expansion  which  gives  to  first  order 


(— ) 
v  3a 


where  (3Go/3a)aQ  is  the  derivative  of  GQ  with  respect  to  an  arbitrary 
variation  in  the  parameters  a  and  evaluated  at  aQ.  We  can  now  write  the 


residual  gradients  from  equation  (37)  as 


-A 

'  3a  ’ 


or  more  explicitly 

-SG  =  ^ +  +  «• *  {k\ se* *  * 


Now  since  there  are  only  five  independent  elements  in  both  Gq  and  5G  we 
have  a  set  of  five  simultaneous  equations  in  five  unknowns  which  can  be 
easily  solved  through  the  use  of  a  standard  computerized  matrix  inversion 
routine,  to  yield  the  corrections  to  the  coil  dipole  parameters,  6<f>r,  50^, 
64,  60m,  and  fig  which  will  eliminate  the  residual  gradients.  The  com¬ 
puter  program  which  performed  these  calculations  called  COILG,  is  des¬ 
cribed  in  Appendix  D 


To  perform  the  coil  alignment  we  first  perform  a  set  of  gradio- 
meter  rotations  in  the  gradient  field  of  the  tower,  find  the  equivalent 
dipole  using  the  computer  programs  GRADB  and  DIPOL  and  align  the  coils 
to  cancel  it.  We  then  perform  a  second  set  of  rotations  and  use  GRADB 
to  compute  the  residual  gradient  matrix,  6G.  The  parameters  specifying 
the  equivalent  dipole  of  the  tower  and  the  residual  gradients  are  then 
used  in  COILG  to  compute  corrections  to  the  dipole  parameters  for  the 
coils  and  the  coil  position  and  currents  are  adjusted  accordingly  to 
cancel  the  residual  gradients.  The  process  can  then  be  repeated  until 
the  gradients  are  reduced  to  a  level  comparable  to  the  inaccuracies  of 
the  gradiometer  measurements;  that  is,  until  the  process  ceases  to  con¬ 
verge. 

Since  the  technique  was  not  developed  until  late  August  we  did  not 
have  an  opportunity  to  test  a  full  sequence  of  rotations  to  complete  the 
iterative  process.  Vie  did  perform  a  sequence  of  three  rotations  consist¬ 
ing  of  a  reference  measurement  with  the  coils  off,  the  initial  measure¬ 
ment  with  the  coils  on  and  another  set  of  rotations  after  adjusting  the 
coils  on  the  basis  of  the  first  two  rotations.  The  results  of  these  measure¬ 
ments  are  shown  in  Table  IV  .  The  required  corrections  are  shown  in  column 
5  where  the  correction  in  the  dipole  position  is  given  by  A0r  and  A<|>r, 
the  correction  in  its  orientation  by  Ad>m  and  A0m  and  Ag  gives  the  required 
change  in  dipole  magnitude.  In  the  first  iteration  the  magnitude  of  the 
residual  gradient  (as  indicated  by  g,  the  magnitude  of  the  equivalent  di¬ 
pole  of  the  residual  gradients)  dropped  from  about  46  NT/m  to  about  23  NT/m. 
Subsequent  iterations  could  undoubtedly  have  improved  on  this  value  but 
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Table  IV. 

Date 

| 

Results  of  calculations  using  COILG  routine  to  determine  cor¬ 
rections  to  coil  dipole  parameters  which  will  nullify  residual 
gradients. 

g  Gradient  elements  Coil  settings  Corrections  to 

(nT/m)  from  GRADB  routine  tL/M  /M  coil  dipole  Param 

9ll/912/913/922/923  (nT/m)  aWab  /a*  / 

(nT/m)  iV4V'V 

“m/A9 

14  Aug  78 

207.9 

-156.43 

0 

reference  rota¬ 

-3.19 

0  • 

tion  to  compote 

5.79 

0 

uncancelled  grad¬ 

344.67 

ients  from  towe i 

-47.9 

‘  | 

15  Aug 

45.7 

13.6 

19.8 

.56° 

4.9 

74.1 

2.16° 

-2.5 

103.0 

8.68°  ;> 

-0.5 

9.52° 

47.5 

26.3  nT/m 

( 

16  Aug 

22.7 

1.0 

-5.1 

-.63° 

-13.2 

-60.5 

1.67°  o 

5.0 

123.0 

-5.27° 

-19.2 

-6.46° 

-3.6 

-6.78  nT/m 

{ 

the  experiment  was  terminated  when  we  began  daytime  collections  of  internal 
wave  data.  The  ultimate  effectiveness  of  the  iteration  with  our  current 
instrument  is  not  clear  from  the  very  few  preliminary  tests  we  performed 
although  the  factor  of  two  reduction  on  the  first  iteration  was  encourag¬ 
ing.  In  view  of  the  problem  outlined  in  Section  VII  regarding  the  appar¬ 
ent  change  in  balance  with  changes  in  the  gradiometer  tilt  angle  and  the 
limitations  on  the  GRADB  calculation  discussed  earlier  in  this  section, 
our  original  goal  of  achieving  a  reduction  of  a  factor  of  100  in  the 
gradient  strength,  g,  is  probably  beyond  the  capability  of  the  instrument. 
However  based  on  our  first  two  measurements  of  residual  gradients  using 
the  C0IL6  approach  we  could  probably  have  achieved  a  somewhat  better  can¬ 
cellation  than  our  best  effort  to  date  which  is  represented  in  Table  III 
by  the  data  recorded  on  5  June.  With  an  improved  instrument,  the  effec¬ 
tive  cancellation  should  be  limited  only  by  the  precision  of  the  gradiom¬ 
eter  measurements  and  our  ability  to  calculate  the  true  residual  grad¬ 
ients,  as  distinct  from  the  effects  of  gradiometer  imbalance. 
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IX.  MODELING  THE  MNO  MEASUREMENT 


In  section  VI  we  discussed  briefly  an  apparent  inconsistency  concern¬ 
ing  the  discrepancy  between  the  MNO  calculated  from  the  parameters  which 
specify  the  equivalent  dipole  of  the  tower  gradients  and  the  MNO  measured 
using  the  gradiometer  noise  spectrum  above  .3  Hz.  We  incorrectly  expected 
that  the  MNO  would  be  that  orientation  described  in  section  V  in  which 

A 

the  dipole  lies  on  the  x-j  axis  of  the  gradiometer  basis  and  the  gradiometer 

A 

pick-up  loops  lie  in  the  plane  defined  by  the  x1  axis  and  the  dipole  moment, 
in,  of  the  equivalent  dipole.  As  outlined  in  detail  in  Podney  and  Gillespie, 
1977  at  this  orientation  the  gradiometer  is  sensitive  only  to  fluctuations 
in  the  polar  angle,  <j>r,  which,  together  with  the  azimuthal  angle,  9r, 
defines  the  position  of  the  equivalent  dipole  with  respect  to  the  gradiom¬ 
eter.  In  other  words,  if  the  gradiometer  output  is  given  by  r  and  the 
dipole  position  by  the  angles,  6r,  4>r ,  its  orientation  by  the  angles, 

0m,  <|>m  and  the  dipole  magnitude  by  g  then  to  first  order 


_ar_  =  _9L_3L_3r_n  _3r_ 

39 „  "  3<f>  30  3g  U’  3<f> 

r  m  m  J  r 


(42) 


We  initially  anticipated  that  these  conditions  would  define  the  orientation  at 
which  the  noise  produced  in  the  gradiometer  by  motions  of  the  instrument 
in  the  large  dc  gradient  field  of  the  tower  would  be  minimized.  However, 
if  3r/3<f>r  is  large  and/or  the  average  fluctuations  in  cjy  are  large  there 
may  be  a  different  orientation  at  which  more  than  one  of  the  derivatives 
of  r  are  nonzero  but  the  resulting  noise  is  still  smaller  than  that  pro¬ 
duced  by  the  fluctuations  in  4>r -  In  view  of  the  discrepancy  between  the 


measured  MNO  and  the  calculated  MNO  we  Investigated  this  question  In  some 
detail  after  the  conclusion  of  the  Internal  wave  measurements  In  August 
of  1978. 


A.  Experimental  Observations 

During  our  measurements  of  the  MNO  during  the  1978  tower  experiments 
we  made  the  following  empirical  observations. 

a.  from  the  systematic,  sets  of  gradlometer  rotations  in  Table  111 
on  31  May  and  5  June  we  found  that  we  could  reduce  the  steady  gradient 
field  by  about  a  factor  of  15  using  the  nested  colls.  This  indicates  that 
the  coils  were  aligned  to  within  about  3°  of  the  true  equivalent  dipole 

of  the  tower.  Based  on  these  measurements  we  expected  the  MNO  to  be  approx i 
mately  at  the  orientation  given  by  the  gimbal  angles: 

$  -  193°  ,  U  -  77°  .  0  -  23°  (43) 

b.  When  measuring  the  MNO  as  described  in  section  VI,  we  initially 
measured  values  of  <j>  and  W  and  then  varied  0  to  find  a  minimum.  Even  on 
this  first  iteration  the  minimum  noise  was  clearly  observed  at  a  value  of 
0  in  the  region  of  7°,  and  not  near  20"  as  expected.  Further  experiments 
confirmed  that  the  noise  level  for  0  near  20"  was  distinctly  greater  than 
for  0  a  7.7°. 

c.  An  overnight  data  collection  with  $  -  197.5°,  U  ■  72.°,  and 

0  ■  23.3°  (the  spectrum  of  which  is  shown  in  Figure  11  as  "coils  off") 
showed  that  at  this  orientation  the  low  frequency  noise  (below  10  mil/.) 
also  was  significantly  greater  than  at  the  measured  MNO. 
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d.  Later  measurements  of  the  MNO  disclosed  an  apparent  degeneracy  in 
the  MNO  defined  by  +  ft  =  270°.  More  specifically,  we  found  that  if  we  held 
0  constant  and  varied  $  and  Q  keeping  <|)  +  n  =  270°,  there  was  no  apparent  in¬ 
crease  in  the  noise  in  the  frequency  range  of  .5  to  10  Hz. 

B.  Model  for  Computing  Gradiometer  Noise 

To  investigate  the  above  observations  in  more  detail  we  attempted  to 
model  the  gradiometer  output  during  MNO  measurements  at  the  tower  in  the  fol¬ 
lowing  manner.  We  first  specify  the  five  parameters  of  the  equivalent  dipole 
of  the  tower  and  compute  the  five  gradient  elements  produced  by  this  dipole 
at  the  gradiometer.  We  then  specify  small  perturbations  in  the  five  dipole 
parameters,  compute  the  corresponding  perturbations  in  the  gradient  elements, 
and  finally  calculate  from  equation  (11)  the  resulting  gradiometer  response 
for  various  gradiometer  orientations.  The  gradiometer  noise  power  is  computed 
as  the  sum  of  the  squares  of  the  perturbations  produced  by  the  individual  di¬ 
pole  parameters  as  written  in  equation  (42). 

The  five  dipole  parameters  are  specified  as  follows: 

<f>r,Gr:  Polar  and  azimuthal  angles  in  the  earth's  basis  specifying  the 
direction  of  the  vector,  r,  from  the  gradiometer  to  the  equiva¬ 
lent  dipole. 

g:  Strength  of  the  equivalent  dipole  given  by  equation  (18). 
x:  Angle  between  the  vector,  r,  and  the  moment,  m,  of  the  equivalent 
dipole. 

A  ->■ 
w:  Angle  between  the  x3  axis  of  the  earth  basis  and  the  vector  rxm. 

The  angles  x  and  <*>  specify  the  dipole  orientation  with  respect  to  the  posi¬ 
tion  vector,  r.  Although  we  could  also  have  defined  the  dipole  orientation 


by  its  polar  and  azimuthal  angles  in  the  earth/s  basis,  the  x.  w  representa- 

j 

tion  has  the  following  advantage.  As  shown  in  section  IV,  equation  (17), 
the  gradients  in  a  reference  frame  fixed  to  a  dipole  are  completely  specified 
by  the  angle  \»  and  the  gradient  strength,  g.  Consequently,  we  can  model 
motion  of  the  gradiometer  with  respect  to  the  tower  by  holding  g  and  \  fixed 
and  varying  #r,  0^,  and  w.  for  a  more  detailed  representation  we  can  select 
perturbations  of  t|>r,  0^,  and  w  characteristic  of  motions  of  the  boom  and  sim¬ 
ultaneously  choose  perturbations  of  \  and  g  characteristic  of  fluctuations 
in  the  tower  dipole  expected  from  ionospheric  effects. 

C.  Comparison  of  Model  with  Experimental  Observations 

The  computer  program  we  use  to  perform  the  calculation  has  been  dubbed 
GNOISE  and  is  described  in  Appendix  E.  In  Table  V  we  show  some  of  the  results 
of  our  modeling  attempts.  In  this  case  we  have  selected  values  of  Ai^,  A0(., 
and  Aw  commensurate  with  the  measured  stability  of  the  cantilevered  boom  of 
about  5  arc  seconds  (as  measured  by  a  biaxial  tiltineter  attached  to  the  grad¬ 
iometer  gitnbal)  and  a  A\  which  we  might  expect  from  fluctuations  of  the  tower 
dipole  driven  by  variations  in  the  earth's  ionosphere.  The  value  of  Ag  was 
selected  to  represent  changes  in  the  effective  dipole  strength  due  to  changes 
in  the  distance  between  the  gradiometer  and  the  tower  dipole.  Specifically, 
the  dipole  parameters  and  perturbations  used  were 


(1>  = 

r 

88° 

A< P  =  .001 

0  = 

110° 

A0  -  . 001 

r 

r 

10  = 

278.7° 

Am  =  . 0006 

X  = 

138.1° 

AX  =  .001 

9  * 

208  nT/in 

Ag/g  =  .0001 
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Table  V.  Output  listing  of  program  GNOISE  which  models  gradlometer  noise  as  a  function  of  gradion- 
eter  orientation.  First  page  lists  dipole  parameters  with  resulting  steady  tower  grad¬ 
ients  and  variation  in  gradients  with  variation  of  dipole  parameters.  Initial  angles 
and  increments  specify  range  of  gradiorreter  orientations  to  be  considered.  Following 
pages  list^the  square  of  graoiometer  output  for  specified  dipole  fluctuations  in  (pT/m)2 
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and  the  corresponding  gradiometer  noise  power  is  printed  out  in  units  of  (pT/m) 
as  a  function  of  all  three  gimba;  angles  <f>,  ft,  and  0. 

The  most  obvious  feature  of  the  calculations  is  that  for  each  value  of 
6  there  is  a  distinct  locus  of  points  described  by  the  constraint  $  +  ft  =  270° 
along  which  the  computed  gradiometer  response  is  a  minimum.  Furthermore,  the 
minimum  in  the  gradiometer  response  as  a  function  of  0  is  at  0  =  9.5°  (not 
at  9  s  23°  as  initially  expected).  This  observation  is  more  clearly  demon¬ 
strated  in  Table  VI  where  we  plot  the  computed  gradiometer  noise  power  as  a 
function  of  all  three  gimbal  angles  over  a  finer  grid.  It  is  clear  that  the 
absolute  minimum  is  at  0Q  =  9.5°,  4>Q  =  235°  and  ftQ  =  35°.  Although  these  values 
for  6  and  ft  differ  significantly  from  our  measured  values  of  212.8°  and 
58.5°  respectively,  nonetheless,  the  measured  values  of  <j>0  and  ftQ  give 

<f>  +  ft  =  271.3°  (45) 

in  excellent  agreement  with  the  computations.  (The  1°  error  in  the  measured 
value  of  <{>0  +  ftQ  is  probably  due  to  a  slight  misalignment  of  the  gradiometer 
with  respect  to  the  gimbal  axes.  This  alignment  was  checked  on  the  tower 
during  the  May,  1978  experiment  but  the  accuracy  with  which  the  alignment  can 
be  performed  is  of  order  1°.)  Also  note  that  the  magnitude  of  the  computed 
gradiometer  noise  power  in  Table  V  shows  that  for  a  given  value  of  0,  the 
gradiometer  noise  power  changes  very  slowly  along  the  diagonal  defined  by 
<t>  +  ft  =  270°. 

In  short,  we  assumed  that  the  tower  gradients  in  Table  III  measured  on 
31  May  approximately  represent  the  equivalent  dipole  of  the  tower  as  given 
by  the  five  parameters  specified  in  section  B  above.  Using  uniform  independent 


Table  VI.  Output  listing  of  program  GNOISE  which  models  gradiometer  noise  as  a  function  of  gradiom- 
eter  orientation.  First  page  lists  dipole  parameters  with  resulting  steady  tower  grad¬ 
ients  and  variation  in  gradients  with  variation  of  dipole  parameters.  Initial  angles 
and  increments  specify  range  of  gradiometer  orientations  to  be  considered.  Following 
pages  list  the  square  of  gradiometer  output  for  specified  dipole  fluctuations  in  (pT/n)* 
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fluctuations  of  the  five  dipole  parameters  the  model  reproduces  virtually  every 
qualitative  feature  of  our  MNO  measurements  on  the  tower. 

We  located  the  MNO  at  0  ®  7.7°,  <p  c  212.8°  and  i)  -  58.5°  which  closely 
corresponds  with  a  reqion  of  reduced  noise  in  the  computational  approach. 
Furthermore,  we  observed  the  degeneracy  in  <fi  +  (i  *  270 °  but  experimental  limi¬ 
tations  imposed  by  the  Inherent  noise  levels  of  both  our  gradiometer  and  the 
spectrum  analyzer  prevented  us  from  finding  the  absolute  minimum  within  the 
$  ♦  fi  degeneracy.  In  further  experiments  with  improved  instrumentation  how¬ 
ever  the  technique  might  provide  important  information  concerning  the  optimal 
gradiometer  alignment  in  more  complicated  magnetic  environments . 
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X.  CONCLUSION 


The  primary  purpose  of  this  work  was  to  develop  noise  suppression  tech¬ 
niques  which  would  allow  us  to  measure  fluctuating  magnetic  gradients  of 
order  1  pT/m  in  the  presence  of  a  200  pT/m  steady  gradient  field.  Gradient 
noise  in  our  geometry  arises  from  magnetization  and  eddy  currents  in  the 
NOSC  tower,  and  small  motions  of  the  gradiometer  in  the  steady  gradient  field. 
From  the  reference  frame  of  the  gradiometer  we  represent  all  of  these  effects 
as  small  fluctuations  in  the  steady  gradient  field  of  the  tower.  We  can 
also  represent  the  gradients  at  any  point  in  space  as  arising  from  a  simple 
magnetic  dipole,  called  the  equivalent  dipole,  and  treat  small  fluctuations 
in  steady  gradients  as  perturbations  in  the  position  and  orientation  of  the 
dipole.  Hence  we  represent  gradient  noise  in  the  gridiometer  as  resulting 
from  small  motions  of  the  equivalent  dipole  of  the  steady  gradient  field. 

The  symmetry  of  this  simple  picture  suggested  that  at  certain  orienta¬ 
tions  the  gradiometer  sensitivity  to  small  dipole  fluctuations  may  be 
reduced.  To  suppress  gradient  noise  then,  we  look  for  a  "minimum  noise 
orientation"  (MNO)  at  which  the  gradiometer  sensitivity  to  dipole  fluctu¬ 
ations  is  minimized.  Experimentally  the  MNO  is  determined  by  measuring  the 
gradiometer  noise  as  a  function  of  orientation  as  the  gradiometer  experiences 
small  motions  in  the  steady  gradient  field.  We  found  that  aligning  the 
gradiometer  to  the  experimentally  determined  MNO  reduced  the  gradiometer 
noise  essentially  to  instrument  noise. 

We  subsequently  modeled  the  gradiometer  noise  using  a  computer  calcu¬ 
lation  in  which  we  computed  separately  the  gradient  fluctuation  from  a  small 
perturbation  in  each  of  the  five  parameters  specifying  the  equivalent  dipole. 
We  then  represented  the  square  of  the  gradiometer  signal  generated  by  an 
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arbitrary  dipole  perturbation  as  the  sum  of  the  squares  of  the  individual 
fluctuations,  and  displayed  this  calculated  gradiometer  "noise"  as  a 
function  of  orientation.  The  results  of  this  calculation  reproduced 
essentially  all  features  of  our  experimental  observations  and  the  calculated 
MNO  agreed  well  with  the  empirically  measured  MNO.  We  conclude  then  that 
our  instrument  responds  as  expected  to  small  gradient  fluctuations,  and 
that  we  can  adequately  represent  small  gradient  fluctuations  as  perturbations 
in  the  position  and  orientation  of  the  equivalent  dipole. 

Although  the  experimental  determination  of  the  MNO  does  not  require 
a  priori  knowledge  of  the  ambient  gradient  field,  to  obtain  a  first  guess 
for  the  MNO  we  use  a  systematic  set  of  gradiometer  rotations  about  the  three 
independent  Euler  angles  to  map  out  the  steady  gradient  field.  From  our 
model  for  the  gradiometer  response  under  rotations  we  deduce  the  position  and 
orientation  of  the  equivalent  dipole  and  use  these  parameters  to  make  an 
initial  estimate  of  the  MNO. 

Another  important  parameter  of  the  gradiometer  is  its  balance.  Small 
differences  in  area  and  orientation  of  pick-up  loops  generate  an  imbalance 
which  can  produce  gradiometer  noise  either  by  temporal  fluctuations  of  a 
spatially  uniform  field  or  small  gradiometer  rotations  in  a  steady  uniform 
field.  For  our  application,  in  which  we  used  the  gradiometer  on  a  very  rigid 
platform,  balance  considerations  did  not  represent  a  limiting  factor, 
although  we  did  investigate  the  balancing  features  of  the  instrument  in  some 
detail . 

In  addition  to  providing  gradient  information,  a  systematic  set  of 
gradiometer  rotations  should  also  provide  knowledge  of  the  gradiometer  balance 
and  the  relative  magnitudes  of  components  of  any  ambient  uniform  magnetic 
field.  We  find  that  values  deduced  from  rotations  of  our  instrument  for  its 
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balance  and  for  components  of  the  earth's  field  are  obviously  erroneous. 

The  failure  of  the  calculations,  which  are  based  on  our  mathematical  model 
for  the  gradiometer  response  under  rotations,  probably  results  from  anamolous 
behavior  of  the  gradiometer  balance. 

The  apparent  balance  of  our  gradiometer  appears  to  change  as  much  as 
a  factor  of  5  as  the  gradiometer  is  tilted  away  from  a  vertical  orientation. 
Within  tilt  angles  less  than  about  20°,  the  effect  is  apparently  non-hysteretic. 
The  origin  of  the  behavior  is  uncertain  but  we  believe  it  may  originate  in 
temperature  fluctuations  in  the  G-10  walls  of  the  dev/ar  or  changes  in  gradio¬ 
meter  orientation  with  respect  to  the  helium  bath.  A  larger  hysteretic 
effect  is  also  observed  for  larger  tilt  angles  which  we  believe  to 
be  caused  by  a  mechanical  motion  of  some  component  of  the  cryogenic  system 
as  the  dewar  is  tilted.  Since  this  effect  occurred  only  at  tilt  angles 
greater  than  about  20°  (which  is  the.  same  constraint  imposed  by  the  dynamic 
range  of  the  instrument  when  in  the  presence  of  the  200  nT/m  gradient 
field  of  the  tower)  the  restriction  did  not  represent  a  serious  limitation 
on  our  experiments. 

The  hysteretic  effect  we  believe  is  probably  characteristic  of  our 
particular  instrument  and  can  be  avoided  with  proper  gradiometer  and 
cryogenic  dewar  design.  The  non-hysteretic  behavior  however  may  be  common 
to  other  superconducting  gradiometer  systems  as  their  balance  approaches 
a  part  in  106  and  if  so  could  represent  a  significant  obstacle  to  the  use 
of  these  devices  on  nonstationary  platforms. 

In  summary,  for  both  large  angle  rotations  and  small  gradient  fluctu¬ 
ations,  our  instrument  responds  to  gradient  fields  as  expected.  However, 
the  gradiometer  response  to  a  uniform  field  through  imbalance  is  distinctly 
peculiar  and  further  work  with  other  instruments  will  be  required  to  resolve 
some  of  the  remaining  questions. 
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APPENDIX  A 


ANAMOLOUS  BEHAVIOR  OF  GRAD IOMETER  BALANCE 

We  originally  believed  that  the  strange  behavior  of  the  apparent  grad- 
iometer  balance  might  be  explained  if  we  were  making  measurements  in  the  pres¬ 
ence  of  a  local  anamolous  gradient  in  the  earth's  field.  We  have  discounted 
that  possibility  however  by  the  following  argument.  From  Table  I  we  note  that 
for  the  data  from  both  Figures  11  and  12,  the  sin  ft  component  changed  by  about 
.25  nT/m  when  we  changed  0.  From  equation  (15)  this  is  equivalent  to  a  large 
change  in  the  sin  ft  coefficient;  specifically  a  change  in  B(9,4>).  If  we  anal¬ 
yze  B ( 0 ,<+> )  for  =  0  as  for  our  data  we  have: 

B(0,0)  =  g^  sin  0  -  g^  cos  0  +  6^  sin  0  -  6-jb^  cos  0  -  6^  ( A- 1 ) 

now  changing  0  from  15°  to  0°  produces 

AB  =  B(15,0)  -  B(0,0)  (A-2) 

or 

AB  =  g^  sin  15°  -  g^(cos  15°-1)  +  <5^b  sin  15°  -  A^Uos  15°-1) 

(A-3) 

If  we  now  assume  6^  and  A^  are  constant  and  note  that  b^  =  0  and  2b£  =  b^ 
for  our  latitude  then 

AB  =  .  26g  1 2  +  .034g^  +  .556^^ 

AB  =  .26g^  +  .034g^  +  (A-4) 

Now  if  there  are  no  gradients  (9-j 2=^1 3=<^^  anc*  ^1  does  not  c^an9e  ^en  B(0,<f>) 
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can  change  by  roughly  50‘t  at  most,  not  a  factor  of  7  as  in  Table  I.  I  f  we 
now  assume  the  observed  AB  of  .25  nT/m  arises  from  gradients  then  we  must  have 
at  least  either  -  1  nT/m  or  g^  =  .75  nT/m,  or  some  combination  of  the 
two. 

First  consider  g^.  If  we  set  12  s  0,  9  =  15°  and  perform  a  $  rotation 
then  from  equation  15  we  have 

r(0,*,O)  »  C{0,$)  ♦  A(9,*)  (A-5) 

where  A( 0 )  is  given  in  equation  (15).  But  this  expression  contains  a  sin¬ 
gle  Fourier  component  of  sin  2<p  of  the  form  1/2  sin  2$  sin  29,  or  for 
9  =  15°,  1/4  g^,,  sin  2$.  The  results  of  a  4>  rotation  with  O  =  15°  and  SJ  =  0 
are  shown  in  Figure  A-l  along  with  the  Fourier  decomposition.  (The  data  were 
taken  at  the  same  time  as  those  shown  in  Figure  12.)  The  value  of  .082  nT/m 
for  the  sin  2$  component  gives 

gj2  *  .33  nT/m  (A-6) 

in  good  agreement  with  values  measured  on  13,  17,  and  19  Jan.  shown  in  Table 
II.  Hence,  the  g^  gradient  element  is  a  factor  of  three  too  small  to  explain 
all  the  effects  shown  in  Figures  11  and  12.  If  we  indeed  take  g^  as  in  equa¬ 
tion  (A-6)  we  would  still  require  g^  3  5  nT/m  which  is  difficult  to  attri¬ 
bute  to  a  geological  anamoly  and  is  more  than  a  factor  of  10  larger  than  the 
largest  g^3  measured  in  Table  II.  Although  the  data  in  Table  II  indicate 
that  longitudinal  gradients  at  La  Posta  may  be  as  large  as  1  nT/m,  we  note 
from  equations  (A-l)  and  (A-3)  that  longitudinal  gradients  cannot  produce  the 
behavior  depicted  in  Figures  11  and  12. 
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Figure  A- 


G=  15° 

n  =  o° 


time 


sin<P:  -.174 
cos?:  .073 
sin  2  <p:  -.082 
cos  2  <f>:  -.093 


.  Chart  recorder  trace  showing  gradiometer  output  during  a  step¬ 
wise  <|>  rotation  with  0  =  15°,  ft  =  0°.  Data  were  taken  at  same 
time  as  data  in  Figure  12.  Numbered  flat  regions  mark  values 
of  4>  at  dwell  points  where  gradiometer  output  was  recorded.  Spikes 
arise  from  eddy  currents  in  mylar  insulation  during  rotation. 
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APPENDIX  B 


GRADB  ROUTINE 


The  GRADB  routine  performs  the  initial  operation  on  data  collected 
during  a  series  of  rotations  of  the  gradiometer.  The  calculation,  which 
is  based  on  the  model  for  the  gradiometer  response  represented  by  equation 
15,  uses  the  gradiometer  data  and  certain  calibration  factors  as  inputs 
and  calculates  as  outputs  the  five  independent  gradient  elements,  g^, 
g-j 2 *  922’  ancl  ^23’  ^ ree  components  in  the  gradiometer  basis  of 

the  imbalance  vector,  3,  and  the  three  components  of  the  earth's  field. 

The  program  can  also  be  operated  in  a  "test  mode"  in  which  the  five 
gradient  elements  and  the  products  6.b.  are  inputs  and  the  routine  calcu- 

*  J 

lates  the  gradiometer  data  for  a  specified  set  of  rotations  then  uses  these 
data  to  recompute  the  gradient  elements  and  the  components  of  3  and  £q. 

This  procedure  allowed  extensive  testing  of  the  routine  for  internal  con¬ 
sistency.  An  added  feature  when  running  GRADB  in  this  mode  is  the  ability 
to  add  random  noise  to  the  computed  gradiometer  data  to  investigate  the 
effects  of  various  noise  levels  in  the  data  on  the  computed  values  of  the 
output  parameters. 

A.  GRADB  Calculation 

The  gradiometer  response,  from  equation  (15),  is  written  in  the  fol¬ 
lowing  form: 

r(o,4>,n)  =  c(e,<j>)  +  A(o,<i>)cos  n  +  B(9,<j>)sin  n  (B-i) 

J 
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where 


C(0.$)  »  CQ(0)  +  Cj  (0)cos  4>  +  C2(0)sin  $  (B-2) 

A(0  *J>)  =  Ap(0)  +  Aj(0)cos  $  +  A2(0)sin  <p  +  A3(0)cos  2$  +  A4(o)sin  2$ 

B ( 0 ,<$• )  =  Bq(0)  +  B^Ojcos  <>  +  B2(0)sin  ^  +  B^(0)cos  2^>  +  B4(o)sin  2$ 

Finally  each  of  the  C^(0),  A^(o),  and  B ^ ( 0 )  can  be  written  in  the  form 
(where  again  C^(0)  has  no  20  components): 

Aj(o)  =  a9  ♦  a|  cos  0  +  A?  sin  0  +  A?  cos  20  +  A^  sin  20  (B-3) 

so  that  we  now  have  59  coefficients  in  the  form  of  AJ. ,  B^ ,  and  C^.  Indeed 

many  of  these  will  be  zero  and  others  will  be  redundant  so  that  ultimately 

there  are  only  a  sufficient  number  to  specify  the  five  gradient  elements 

and  the  nine  products  S.B.  where  i , j  =  1,  2,  3. 

3  3 

A  sample  of  the  gradiometer  data  is  given  in  Table  B-J  where  each 
block  of  data  represents  the  specified  value  of  theta,  each  row  within  a 
block  represents  a  specific  value  of  phi  from  60°  to  360°,  and  each  column 
corresponds  to  one  value  of  Q  from  90°  to  360°.  The  calculation  then  pro¬ 
ceeds  as  follows.  The  Fourier  decomposition  of  equation  (B-l)  on  each  row 
produces  a  value  of  C(0,<?),  A( 0 )  and  B(0,<t>)  corresponding  to  each  value 
of  0  and  <p.  A  second  Fourier  decomposition  in  41  on  these  sets  of  numbers 
produces  tht  coefficients  A^(0),  B^(0),  and  C^(0).  These  15  sets  are  or¬ 
ganized  according  to  0  as  shown  in  Table  B-II,  and  coefficients  A'? ,  B'j , 
and  are  then  calculated  by  a  third  Fourier  decomposition  in  0.  The 
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Table  B-I.  Sample  of  gradiometer  data  collected  during  a  systematic  set  of 
gradiometer  rotations.  Data  are  grouped  according  to  values  of 
0  as  indicated.  Within  a  block,  each  row  corresponds  to  one 
value  of  <j»  from  60°  to  360°  and  each  column  to  one  value  of  ft 
from  90°  to  360°.  Data  are  recorded  in  units  of  volts  directly 
from  gradiometer  output. 
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Table  B- 1 1 .  Sample  of  Fourier  coefficients  from  Fourier  decompositions  in 

SI  and  <j>  from  data  in  Table  B- 1 .  Correspondence  of  coefficients 
in  text  with  those  shown  here  is  as  follows: 

CQ(e)  =  AOAO,  C]  =  A0A1 ,  C?  =  A0B1 

Aq  ■  A1A0,  Aj  =  A1A1 ,  A^  =  A1B1 ,  A3  =  A1A2,  A4  =  A1B2 

Bq  =  B1A0,  Bj  =  B1A1 ,  B;)  =  B1B1,  B3  =  B1A2,  B4  =  B1B2 
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gradient  elements  and  products  6.B.  are  then  computed  and  printed  out. 

1  J 

A  flowchart  of  the  normal  mode  processing  is  shown  in  Figure  B-l.  Input 
values  for  the  gradiometer  calibration  constant  and  scale  factor,  the  earth's 
field  magnitude,  the  parameters  specifying  the  exact  gradiometer  rotations, 
and  the  gradiometer  data  are  obtained  from  a  disc  file,  $GRAD#,  where 
the  value  of  #  is  specified  by  the  operator. 

B.  Test  Mode 

When  running  the  program  in  the  test  mode,  which  is  activated  by 
an  operator  option  at  program  initiation,  the  program  requests  the  five 
gradient  elements,  the  products  S.b.,  and  the  desired  noise  level  as 

•  vl 

inputs,  then  obtains  the  gradiometer  constants,  earth's  field  magnitude 
and  desired  rotations  from  file  $BLNC0.  The  program  then  computes  the 
expected  gradiometer  response  from  equation  (15)  at  each  value  of  0,  <p,  and  ft, 
adds  the  specified  level  of  random  noise  and  stores  the  data  in  file  $BLNC0. 
The  test  mode  then  performs  exactly  the  same  processing  on  the  artificial 
data  now  residing  in  file  $BLNC0  as  does  the  normal  mode  processing  on 
the  data  in  a  file  $GRAD#.  A  flowchart  of  the  test  mode  operation  of 
GRAD8  is  shown  in  Figure  B-2.  A  listing  of  the  complete  routine  is  in¬ 
cluded  in  Appendix  G. 

We  performed  several  test  runs  of  the  GRADB  routine  where 
we  specified  the  dipole  parameters  along  with  a  random  noise  contribution 
to  the  gradiometer  data.  Even  for  a  random  noise  contribution  which 
represents  fluctuations  of  order  10  percent  in  the  original  data  the  GRADB 
program  produces  the  correct  values  for  the  gradient  elements  and  the 
imbalance  vector  to  within  one  percent. 


C .  GRADB  Inputs 

1.  Operator  inputs  -  normal  mode: 

a.  file  name  $GRAD#  containing  data  to  process  where  0  <  #  <  9. 

2.  Operator  inputs  -  test  mode: 

a.  dc  offset  level  of  gradiometer 

b.  5  gradient  elements  specifying  gradient  field 

c.  9  products  of  components  of  imbalance  vector  3  and  earth's 
field, 

d.  noise  magnitude  to  be  applied  to  calculated  gradiometer  response 

3.  File  inputs  -  both  modes: 

a.  gradiometer  calibration  constant  [v/(nT/m)l 

b.  gradiometer  scale  factor 

c.  magnitude  of  earth's  field  in  nT 

d.  increments  in  0,  <J>,  and  0 

e.  maximum  value  of  0  for  gradiometer  rotations 

f.  gradiometer  data  array  from  specified  file 

D .  GRADB  Outputs 

1.  Normal  Mode: 

a.  dc  offset  level  of  gradiometer 

b.  5  gradient  elements  specifying  gradient  field 

c.  9  products  of  components  of  imbalance  vector,  &,  and  earth's 
field,  Ij. 

d.  3  components  of  3 

e.  3  components  of  t>0 

f.  Fourier  components  of  ft  and  <j>  decompositions  (as  in  Table  B-Il) 

2.  Test  Mode: 

a.  The  test  mode  produces  all  of  the  line  printer  outputs  speci¬ 
fied  above  and  in  addition  generates  artificial  gradiometer 
data  with  the  specified  noise  level  in  file  0BLNCO. 
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APPENDIX  C 


DIPOL  ROUTINE 


The  DIPOL  routine  computes  the  equivalent  dipole  of  the  five  grad¬ 
ient  elements,  g^,  g12»  0-)3»  ^22’  and  923’  calculated  by  the  GRADB  routine 
(see  Appendix  A).  Using  the  gradient  matrix,  G,  specified  by  the  five 
g..'s,  the  routine  uses  a  matrix  inversion  subroutine  to  compute  the  eigen- 
values  of  G  and  the  corresponding  eigenvectors  which  specify  the  principle 
axes  of  the  gradient  field.  From  these  data  the  DIPOL  routine  then  cal¬ 
culates  the  location  and  orientation  of  each  of  the  four  equiva¬ 
lent  dipoles  and  a  corresponding  minimum  noise  orientation. 


A.  DIPOL  Calculation 

Given  a  gradient  matrix,  G,  in  the  earth's  basis  we  wish  to  find 
principle  axes  for  which  the  matrix  is  diagonal;  namely, 

A  =  RpGRp  (C-l) 

where  the  matrix  Rp  specifies  the  rotation  which  brings  the  earth's  basis 
into  coincidence  with  the  principle  axes.  The  matrix  A  has  diagonal  ele¬ 
ments  A.j ,  and  A^  that  are  eigenvalues  of  the  matrix  G  and  are  ordered 

so  that  A,  >  A0  >  A0.  Columns  of  the  matrix  R„  are  the  corres- 
1—2—3  p 

ponding  eigenvectors,  which  specify  principle  axes  of  the  gradient  field 
represented  by  G.  This  calculation  is  performed  by  subroutine  OACV  which 
was  adapted  from  an  existing  routine  in  the  FORTRAN  library  of  the  Lawr¬ 
ence  Berkeley  Lab's  CDC  6600.  The  DIPOL  routine  then  calculates  the  position, 
orientation,  and  strength  of  the  four  possible  equivalent  dipoles  and  the 


Euler  angles  of  the  girnbal  which  will  bring  the  gradiometer  axes  into 
coincidence  with  the  principle  axes.  In  addition,  the  MNO  is  calculated 
which  corresponds  to  the  gradiometer  orientation  in  which  the  equivalent 
dipoles  lies  on  the  xj  axis  of  the  gradiometer  basis  and  the  gradiometer 
pick-up  loops  lie  in  the  plane  defined  by  xj  and  the  dipole  moment,  m, 
of  the  equivalent  dipole.  A  flowchart  of  this  processing  sequence  is 
shown  in  Figure  C-l. 

B.  Test  Mode 

In  the  test  mode  the  DIPOl  routine  accepts  as  operator  inputs  the 
position,  orientation,  and  strength  of  a  dipole.  The  gradients  from  the 
dipole  at  the  site  of  the  gradiometer  are  calculated  and  the  five  grad¬ 
ient  elements  are  then  used  by  the  routine  to  perform  exactly  the  same 
calculation  as  performed  in  normal  mode  processing.  The  test  mode  thus 
allowed  us  to  check  the  computer  code  for  internal  consistency  and  study 
the  variations  in  the  gradient  elements  induced  by  small  fluctuations 
in  the  position,  orientation,  and  strength  of  a  dipole  moment  near  the 
gradiometer.  A  flowchart  of  the  test  mode  processing  is  shown  in  Figure 
C-2.  A  program  listing  for  the  DIPOL  routine  is  given  in  Appendix  G. 

C.  DIPOl.  Inputs 

1.  Nonna 1  mode: 

a.  five  gradient  elements 

2.  Test  mode: 

a.  polar  and  azimuthal  angles  specifying  dipole  position 


C-? 


b.  polar  and  azimuthal  angles  specifying  dipole  orientation 

c.  strength  of  dipole  and  distance  from  gradiometer 
PI POL  Outputs 

1.  Normal  Mode: 

a.  prints  out  input  values  of  g . ■ 

'  J 

b.  eigenvalues  and  eigenvectors  of  matrix  G 

c.  polar  and  azimuthal  angles  specifying  orientation  of  each  eigenvector 

d.  Euler  angles  giving  principle  axes 

e.  location  and  orientation  of  each  of  four  possible  equivalent 
dipoles 

f.  MNO  for  each  equivalent  dipole 

2.  Test  mode: 

a.  prints  out  input  values  for  specified  dipole 

b.  five  gradient  elements  g..  produced  by  the  specified  dipole 

*  J 

c.  prints  out  all  other  outputs  listed  under  normal  mode  pro- 
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APPENDIX  D 


CO I LG  ROUTINE 

The  COILG  routine  was  designed  to  allow  a  precise  alignment  of  the 
nested  electrical  coils  in  order  to  cancel  the  steady  gradient  field  of 
the  tower  using  the  equivalent  dipole  concept.  Using  a  systematic  set 
of  rotations  of  the  gradiometer  to  map  out  the  uncancelled  steady  gradient 
field,  the  GRADB  and  DIPOL  routines  outlined  in  Appendices  Band  C  pro¬ 
vide  an  estimate  of  the  position,  orientation,  and  strength  of  the  equiv¬ 
alent  dipole  of  the  steady  gradients.  If  the  coils  are  now  aligned  and 
energized  to  cancel  the  steady  gradient  field  a  second  set  of  gradiometer 
rotations  will  map  out  the  residual  gradients  from  which  COILG  will  esti¬ 
mate  the  corrections  to  the  coil  position,  orientation  and  strength  re¬ 
quired  to  nullify  the  residual  gradients. 

A.  GOILG  Calculation 

The  GOILG  routine  accepts  operator  inputs  specifying  the  position, 
orientation,  and  strength  of  the  coil  dipole  and  the  residual  gradient 
elements  calculated  by  the  GRADB  routine  from  the  second  set  of  gradiometer 
rotations.  From  the  parameters  specifying  the  coil  dipole  the  program 
initially  calculates  the  gradient  matrix,  Gq,  at  the  position  of  the 
gradiometer.  From  equation  (41 )  in  section  VIII  we  can  compute  the  required 
dipole  corrections  from  a  knowledge  of  the  residual  gradients,  6G,  and 
the  derivatives  of  Gq  with  respect  to  variations  in  the  various  dipole 
parameters.  To  calculate  the  derivatives  of  GQ  we  increment  one  of  the 
dipole  parameters,  say  <b  ,  by  some  amount  A<J>r  and  compute  a  new  gradient 
matrix,  G^,  for  4>r  =  +  Wr-  The  derivative  of  G0  with  respect  to  <j>r 

is  then  given  approximately  by 


«G  G ' -G 
3  0  _  0  0 


34> 


A<t> 


(0-1) 


and  similarly  for  the  other  four  parameters  of  the  dipole.  We  now  write 
equation  (41): 
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or  where  G  and  its  derivatives  are  the  matrices: 
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and  the  residual  gradient  elements  are: 


6G  = 


Si 

«9,2 

6913 

«922 

\S923  / 


(0-4) 


Since  6G  and  the  derivatives  of  GQ  are  known,  the  system  of  equations  is 

easily  solved  using  a  matrix  inversion  routine  to  provide  values  for 

6*  ,  66  ,  6<b  ,  60  and  6o.  A  flowchart  for  the  COILG  processing  sequence  and 
r  m  m  * 

a  listing  of  the  fortran  program  are  included  in  Figure  0-1  and  Appendix  G 
respectively. 


B.  COILG  Inputs 


a.  Five  parameters  specifying  the  coil  dipole  moment 

b.  Five  residual  gradients  calculated  by  the  GRADB  routine  from  the  second 
set  of  gradiometer  rotations 

c.  Coil  settings  in  nT/m. 

d.  Perturbations  to  be  applied  to  dipole  parameters  for  calculating 
derivatives  of  GQ(A<t>r,  AO^,  A<t>m>  A6m,  Ag) 

C.  COILG  Outputs 

a.  Variation  of  Gq  with  specified  perturbations  in  dipole  parameters 

b.  Derivatives  of  G 

o 

c.  Corrections  in  dipole  parameters  required  to  nullify  residual 
gradients. 

d.  New  current  settings  for  nested  coils. 


Al’l’l  NDIX  l 


liNOlSl  ROUUNl 


The  GNOISE  routine  calculates  the  gradiometer  response  to  small  fluc¬ 
tuations  in  each  of  the  five  parameters  specifying  an  equivalent  dipole. 
The  total  gradiometer  noise  is  then  computed  as  the  sum  of  the  squares 
of  the  individual  fluctuations  over  a  range  of  gradiometer  orientations 
specified  by  the  operator. 

A.  GNOISE  Calculation 


The  GNOISE  routine  accepts  operator  inputs  specifying  the  position, 
orientation,  and  strength  of  the  equivalent  dipole  and  the  amount  by  which 
each  parameter  is  to  be  perturbed.  The  program  also  requests  the  range 
of  gimbal  angles  specifying  the  orientations  over  which  the  noise  is  to 
be  computed. 

First  the  gradient  elements  at  the  position  of  the  gradiometer  are 
calculated  from  the  five  dipole  parameters  specified  by  the  operator. 

We  write  the  five  independent  gradient  elements  as  a  matrix  G  given  by 


The  program  now  calculates,  <SG,  the  variation  in  G  with  respect  to  each 
of  the  dipole  parameters  as: 


r-i 


(E-2) 


6G  =  ( 6G[4>r ] ,  6G[0r],  6G[<J)m],  AGCoJ,  6G[g]) 


where 


6G[a]  = 


\ 


Aa 


(E-3) 


Now  using  equation  (11)  the  gradiomete^  response  to  a  fluctuation  in  each 
of  the  dipole  parameters  is  computed  as: 

Ara(0,<$>,ft)  =  [rle,<t>,Ji,G)  -  r(0,(|>,ft,G+SG[a])]/Aa  (E-4) 

and  the  computed  gradiometer  noise,  Ng,  is  represented  as 

N  (0,*,n)  =  l  [Ara(0,<|>,ft)]2  (E-5) 

y  a 

The  noise,  Ng,  is  then  printed  in  matrix  format  for  each  0,  <J> ,  and  ft 
over  some  range  of  each  of  the  angles.  A  flowchart  of  the  GNOISE  proces¬ 
sing  sequence  is  shown  in  Figure  E-l  and  a  listing  of  the  fortran  program 
appears  in  Appendix  G. 

B.  GNOISE  Inputs 

a.  Five  parameters  specifying  dipole  moment 

b.  Perturbations  to  be  applied  to  each  of  the  dipole  parameters. 


I.. 


«~jj 
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c.  Initial  gimbal  angles  of  gradiometer. 

c.  Increments  in  gimbal  angles  at  which  the  noise  is  to  be  computed. 


C.  GNOISE  Outputs 

a.  Prints  out  input  parameters. 

b.  Gradients  produced  by  the  specified  dipole. 

c.  Variation  and  derivatives  of  gradients  with  respect  to  varia¬ 
tions  in  each  of  the  dipole  parameters. 

d.  Initial  gimbal  angles  and  desired  increments. 

e.  Gradiometer  noise  for  each  combination  of  0,  <j>,  £2  over  the  range 
desired. 
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APPENDIX  F 


GTILT  ROUTINE 

The  GTILT  routine  is  a  small  utility  program  used  primarily  during 
our  attempts  to  align  the  gradiometer  with  the  ambient  gradient  field. 
GTILT  uses  as  inputs  the  initial  orientation  of  the  gradiometer  as  speci¬ 
fied  by  the  three  Euler  angles  0 ,  ,  and  SI,  an  arbitrary  rotation  axis 

specified  with  respect  to  either  the  earth  or  the  gradiometer  basis,  and 
a  rotation  angle,  It  then  calculates  the  new  Euler  angles  for  the 
gradiometer  which  will  produce  a  rotation  of  the  gradiometer  through  the 
angle,  ip ,  about  the  specified  axis  of  rotation.  Figure  F-l  shows  a  flow¬ 
chart  of  the  program  and  a  listing  is  given  in  Appendix  G. 

A .  GTILT  Inputs 

1.  Euler  angles  specifying  initial  gradiometer  orientation. 

2.  Rotation  axis  specified  with  respect  to  either  the  earth  or  grad¬ 
iometer  basis. 

3.  Rotation  angle,  4>. 

B .  GT ILT  Outputs 

1.  Prints  out  input  values. 

2.  New  Euler  angles  which  produce  the  desired  rotation. 


F-l 


GTILT  FLOWCHART 


Figure  F-l 
F-2 
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PROGRAM  LISTINGS 
GRADB  ROUTINE 
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